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Abstract 

We present an elementary introduction to the problem of quark-hadron duality 
and its practical limitations, in particular as it concerns local duality violation in 
inclusive B meson decays. We show that the accurate definition of duality violation 
elaborated over the recent years allows one to derive informative constraints on vio- 
lations of local duality. The magnitude of duality violation is particularly restricted 
in the total semileptonic widths. This explains its strong suppression in concrete 
dynamical estimates. We analyze the origin of the suppression factors in a model- 
independent setting, including a fresh perspective on the Small Velocity expansion. 
A new potentially significant mechanism for violation of local duality in T S \(B) is 
analyzed. Yet we conclude that the amount of duality violation in T S \(B) must 
be safely below the half percent level, with realistic estimates being actually much 
smaller. Violation of local duality in Y S \{B) is thus far below the level relevant to 
phenomenology. We also present a cautionary note on the B — > D* decay amplitude 
at zero recoil and show that it is much more vulnerable to violations of quark-hadron 
duality than T S \(B). A critical review of some recent literature is given. We point 
out that the presently limiting factor in genuinely model-independent extraction of 
V c b is the precise value of the short- dist ance charm quark mass. We suggest a di- 
rect and precise experimental check of local quark-hadron duality in semileptonic 
B — ► X c Iv decays. 
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1 The problem 

In QCD one evaluates transition rates, distributions etc. in terms of quarks and 
gluons, for which the fundamental interactions are specified. Yet the asymptotic 
and thus observable states are hadrons. Some notion of quark-hadron duality - or 
duality for short - has underlied many applications of the quark model from the 
early days on. It is based on the idea that a quark level calculation should at least 
approximate hadronic rates. Typically it was invoked in a rather vaguely formulated 
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way; for example it was not stated how accurate such an approximation would be. 
A more concrete phenomenological formulation was given in the special case of deep 
inelastic lepton-nucleon scattering to extend the onset of scaling to lower momentum 
transfers |J. 

Nowadays there is little doubt that QCD is the true theory of strong interactions 
predicting - in principle - all properties of hadrons with unlimited accuracy. The 
assumption of duality in its most general version amounts to no more than this 
conviction; i.e., the true hadronic observables coincide with what one ultimately 
obtains in the quark-gluon language provided all possible sources of corrections to 
the parton picture stemming from QCD itself are properly accounted for. The 
question of duality thus has shifted to a different level: it amounts to assessing the 
potential magnitude of those QCD contributions to the observables of interest that 
have not been included. 

The practical validity of duality thus depends on the theoretical tools available 
for treating QCD dynamics. The inception of QCD as a theory of strong interac- 
tions was intimately related with realizing its asymptotically free nature. Corre- 
spondingly, at first the quantitative treatment of strong interactions relied almost 
exclusively on perturbative expansion based on the smallness of the effective gauge 
coupling between quarks and gluons in processes at small space-time intervals. It 
was realized and checked in experiment that the parton ansatz improved by calcula- 
ble perturbative corrections yields a good approximation for a wealth of high-energy 
Euclidean observables representing true short-distance physics. 

Generic QCD effects for real processes in Minkowski space are quite different 
and more complicated. Even at arbitrary high energies they can acquire diver- 
gent corrections already in the perturbative expansion, effects governed not by the 
small running coupling evaluated at the scale of overall energy, etc. A number of 
"infrared-stable" [fj] Minkowskian observables were identified which are free from 
these complications, thus being candidates for applying the concept of duality in 
practice. Those included sufficiently inclusive processes which combined different 
quark-gluon channels in a certain way. Likewise, quark-hadron duality implied to 
be applicable if a sufficiënt number of hadronic channels were included. Yet, once 
again, it remained somewhat indefinite what shall constitute a sufficiënt number of 
channels. 

One of the most inclusive and a priori infrared-stable Minkowskian observable 
amenable to perturbation theory is the total cross section of e + e~ annihilation into 
hadrons.Q The concept of duality there was first addressed theoretically in Refs. || || 
and the more specific notion of local quark-hadron duality was introduced. 

The parton ansatz yields an energy-independent ratio 



1 We assume here annihilation to light-flavor hadrons only. Production of heavy quarks has 
some peculiarities making discussion of duality less transparent. 
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equal to the number of colors N c (hereafter we omit simple factors reflecting quark 
charges, number of navors, etc). The perturbative corrections slightly modify this: 



Yet those are governed by the small running coupling a s (s) and thus are not sig- 
nificant. Moreover, the parton cross section remains smooth and monotonic in 
perturbation theory. 

The experimental cross section, however, exhibits manifest resonance structure 
up to relatively high energies. In this domain duality between the QCD-inferred 
cross section and the observed one at a given fixed energy looked problematic. It 
was suggested that the equality between the two is restored if averaging over an 
energy interval - or 'smearing' - is applied. The problem of how to compare the 
QCD-based and the actual hadronic probabilities point-to-point in energy (or other 
kinematic variables) is referred to as local quark-hadron duality. 

Resonance physics, at least in the light quark sector is tightly related to confining 
properties of QCD. Since the latter has nonperturbative origin, analyzing dynamic 
aspects of local duality requires control over nonperturbative effects in a consistent 
QCD-based framework. 

No real progress beyond the qualitative stage outlined in Refs. |3], [| occurred for 
a long time. In most hard processes the problem could be evaded by going to higher 
energy scales where duality violations are greatly reduced. In beauty decays this 
option does not exist; yet data were of less than sterling quality and therefore did 
not create pressure for a more precise theoretical treatment. On the theoretical side 
there was an unsurpassed road block: as long as one has very limited control over 
nonperturbative effects, there is little meaningful that can be said about duality 
violations. 

The general ideas of the Operator Product Expansion (OPE) have been applied 
to quantify nonperturbative effects in a number of important QCD processes since 
the early 1980's Subsequently data on beauty decays became quite precise as 
well. It was also realized that heavy flavor physics had great potential for revealing 
the presence of New Physics indirectly; such a goal places a premium on accuracy. 
At the same time heavy quark theory was developed allowing to deal with many 
nonperturbative aspects in heavy flavor decays in a well-defined way Those devel- 
opments refocussed attention on the need to deal with duality and its limitations in 
a quant it at ive way. 

While we have no complete theory of duality violation yet, significant progress 
towards that goal has been made, specifically for those Minkowskian observables that 
are described by the OPE. Violation of local duality there is related to the asymptotic 
nature of the power expansion employed by the "practical version" of the OPE P, [71 . 
Such series are not sensitive to contributions which asymptotically behave, say, like 
exp(— constv/QV/ihadr)- While they are exponentially suppressed in the Euclidean 
domain, they turn into pure oscillations upon continuation to Minkowski kinematics, 
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as in the above example with Q 2 = — s = —E 2 m , and remain only power-suppressed 
in absolute magnitude due to various 'preexponential' factors. This is the way how 
the general question of quark-hadron duality inherent in any QCD expansion gives 
rise to local duality violations specific to Minkowskian kinematics. 

At the same time, the OPE imposes essential constraints on possible (local) 
duality violations, which are often missed in the literature. As a result, many effects 
alleged to generate violations of local duality, especially in the total decay rates of 
heavy navors would actually signify violations of the OPE. Thus they have little 
relevance to specifics of the Minkowskian domain, but rather contradict general 
principles of QCD common to any type of power expansion both in Minkowskian 
and Euclidean kinematics. 

Some clear concepts for the physical origin behind duality violations as well as 
their mathematical portals have been identified in recent years. In particular, it was 
shown H that violations of local duality is a very general phenomenon not directly 
associated with confinement or the physics of narrow resonances. A review of our 
present theoretical understanding of local quark-hadron duality and its violation can 
be found in the recent dedicated publication by Shifman |J. 

We have to note at this point that our terminology slightly differs from the one 
adopted in Ref. || dedicated specifically to local duality. In particular, the intrinsic 
limitations on the accuracy of the 'practical' OPE for truly Euclidean observables 
at finite mass scales are not referred to there as related to the general issue of quark- 
hadron duality. As explained in a clear way in the review 0, they are typically self- 
manifest in the OPE itself, and thus can be viewed as expected "natural" or "usual" 
uncertainty. Local duality violation for Minkowskian observables is a phenomenon 
over and above this. For this reason M. Shifman does not put such an emphasis on 
distinguishing the terms local duality violations and duality violation in general, and 
rather often uses the both on the parallel footing. Consequently, in the terminology 
consistently adopted in Ref. || there is by defmition no duality violation in Euclidean 
observables, but only 'natural' uncertainties of the asymptotic power expansions. On 
the other hand, we refer to the latter as limitations of generic quark-hadron duality, 
with local duality violation being its very specific Minkowskian aspect. Since the 
focus of all the discussions lies in local duality where our terminologies are identical, 
the difference should not lead to any confusion in the practical aspects. 

Quantitative tests of local duality violations are ultimately provided by data. Yet 
this is not completely straightforward. For there are practical uncertainties in OPE 
predictions due to uncertainties in the input parameters, which are intrinsically 
unrelated to duality. E.g., variations in the values of the strong coupling, quark 
masses, the leading nonperturbative condensates etc. often lead to more significant 
uncertainties. The number of clean testgrounds for local duality is thus limited in 
practice. 

A striking demonstration of the confidence the HEP community has in the 
asymptotic validity of duality was provided by the discussion of the heavy flavor 
widths of the Z° resonance. There was an about 2% difference in the predicted and 
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the observed Z° decay width into beauty which led to a lively debate on how signifi- 
cant that was vis-a-vie the experimental error. No concern was expressed about the 
fact that the Z° width into beauty hadrons was calculated for quarks, yet measured 
for hadrons. Moreover, the strong coupling a s (Mz) is routinely extracted from the 
perturbatively computed hadronic Z° width with a stated theoretical uncertainty 
0.003 which translates into a confidence in i\ad(-Z°) of about 10 -3 . 

This confidence derives from the fact that duality actually represents a very 
natural concept M. It is based on the picture that OPE-treatable processes with 
hadrons evolve in two steps. At first hard dynamics proceed in the femto universe 
characterized by large scales like mq or momentum transfers \[Q 1 . Subsequently the 
quarks (and gluons) transmogrify themselves into hadrons; since this transformation 
is driven by soft dynamics, this second step is characterized by much larger distance 
scales. Let us in particular consider the decay of a heavy quark into light flavors. 
The typical time scale for the first step is provided by 1/rriQ, for the second one 
by 1/Aqcd m the restframe of a final state quark and thus by rriq/ Aq CD in the 
original rest frame. When the second step occurs, the quarks originally present are 
far removed from each other. One then expects that the second step will determine 
the composition of the final state, but not gross characteristics like total rates, the 
directions of energetic jets etc, since those are well established by that time. 

This can be illustrated in a simple quantum mechanical example. Consider the 
weak decay of a heavy quark Q bound to an antiquark q by a potential V(R) centered 
on Q. The overall decay width in the 1/rriQ expansion is determined by the local 
properties of the potential, namely by V(R) for R ^ ö{l/mQ). Yet the spectrum of 
the final states and other more detailed properties depend - even at large ulq - on 
the details of the potential, like its behavior at finite (or even infinite) distances 
Consequently, certain effect of the long distance properties of the potential must be 
observed at finite rriQ in the total decay width as well - these affect local duality at a 
given rriQ. Applied to actual heavy quarks, the real argument is about how good an 
approximation duality represents at lower scales, in particular at the beauty mass 
scale, and whether it ceases to be of a numerical value at the charm scale or not. 

This Vademecum is organized as follows. After illustrating the problem of local 
duality with the example of the r width, we give our Executive Summary on local 
duality in heavy flavor decays in Sect. 3. Those results are explained in subsequent 
sections. After introducing the theoretical arsenal in Sect. 4, we briefly describe two 
implementations of nonperturbative dynamics in Sect. 5. This framework is applied 
to a dedicated discussion of the total semileptonic width of B mesons in Sect. 6 
complemented with selected comments on the literature in Sect. 7. In Sect. 8 we 
present remarks on the role of duality in describing B — > D* iv at zero recoil before 

2 For example, if V(R) — ► co for R — > co, the q'q spectrum in [Qq] — > Ivq'q will be discrete 
consisting of a series of single narrow mesons; if on the other hand the potential saturates for 
finite R - modcling the possibility for the potential string between q' and q to break producing 
additional quark pairs as it happens for real QCD - the spectrum will be continuous with many 
different hadronic configurations. 
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listing conclusions in Sect. 9 and offering final observations in Sect. 10. 



2 The scope 

While the practical aspects of general quark-hadron duality are involved in any 
dynamical computation one undertakes in QCD - even the more specific aspect of 
local duality has to be faced in a wide range of problems - we limit the discussion 
here mainly to weak decays of heavy navors, with inclusive widths there being one 
of the best studied subjects in this respect. Moreover, we focus on semileptonic 
decay widths from which \V c b\ and \V u b\ are primarily extracted and where duality 
violations are most constrained. To give a sense of the peculiarities of local quark- 
hadron duality we compare this with the zero-recoil B — » D* transition amplitude 
which is closer to a Euclidean-type observable. As expected on general grounds if 
rrib ~ m c were to hold B — > D* would be less vulnerable to duality violation than 
T S \(B). Yet since the energy release in b — > clv is much larger than the charm 
quark mass, the limitations of local quark-hadron duality in r s i(5) are suppressed 
compared to the intrinsic uncertainties of the power expansion for the B — > D* 
amplitude. The hadronic width of r leptons is considered to illustrate the numerical 
aspects of local duality in a simple setting. 



2.1 The hadronic r decay width 

The observable 

D _ r(r~ — > v T + hadrons) 

K r = ^TT-I — T {O) 

is the closest analogy to the integrated - or smeared - cross section of e + e~ anni- 
hilation in the interval of energy up to m T ~ 1 . 77 Ge V. More precisely, up to small 
electroweak radiative corrections it is expressed in terms of the spectral densities py 
and pa in the vector and axial-vector channels: 

R ~ [ M? dS (l S Vfl 4- 23 )\o (* + o fsïl - /o(M ' } 3 /2(M ' } +2 /3(M ' } 

Rt -j ~ïpX~k) [ + m) [pv{)+pAi )] ~^r~ 3 ^r + ^r 

(4) 

(here quark masses are neglected), where the moments /„ are defined as 

I n (M) = / dss n [ Pv (s)+p A (s)} . (5) 

In spite of being an averaged spectral density, R T is affected by physics underly- 
ing local duality, which can be illustrated by the following simple argument. Let us 
consider the chiral limit m u ^, s = and suppose the perturbative effects can be ne- 
glected completely, including all the anomalous dimensions of vacuüm condensates. 
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Eq. (||) would then seem to state - without any reference to strong dynamics - that 
R T obtained in the "practical" OPE is at most a fourth-order polynomial in 1/M 2 : 

p A A2 As Ai (r\ 

R T = A » + MÏ + Ml + M-y (6) 

T T T 

at least above a certain mass scale M . However, since hadronic thresholds in anni- 
hilation exist at arbitrary high energy, this certainly cannot be true. The component 
of R T not contained in Eq.flC]) represents a violation of local duality. While naively 
almost arbitrary, in reality it obeys severe constraints stemming from the OPE. An 
explicit model was analyzed in Ref. which manifestly exhibited these features. 

A critical look at the available experimental data on pv(s) and Pa(s) suggests 
that violation of local duality a priori could be quite sizeable. This is not surprising 
since nearly half of the decay probability is due to hadronic states with invariant 
masses not exceeding 1 GeV, a resonance rather than asymptotic regime. An at- 
tempt is often made to extract a precise value of the strong coupling a s (M T ) from 
R T , neglecting the potential violation of local duality. This is hardly justified the- 
oretically JTT] . However, the numerical estimates of duality violation in R T yield 



rather small values, at a few percent level flT(J , which may sound surprising in view 
of the relatively low momentum scale characteristic to r decays. This is indirectly 
confirmed by experiment. 

The OPE yields the following large-M T expansion of R T : 

R T j-, A / s C2 C3 C4 

-=i?o + A pcrt K) + — + — + — , (7) 

where R — 1 — 4(m 2 +m^ + sin 2 9c{jn 2 s — m^))/M 2 + ... is the parton expression, 
A pert (a s ) represents the perturbative series, and the last terms emerge from the 
nonperturbative condensates. The and 1/M® terms appear only due to the 

presence of perturbative corrections or/and contain powers of the light quark masses. 
As a result, the OPE power corrections turn out to be strongly suppressed, at a 
percent level. 

The perturbative factor A pert (a: s ) = 1 + a s /n + ... has been well studied theo- 
retically. Assuming the canonical values of the condensates, Eq. (0) reproduces the 
experimental value of R T at a s (M T ) ~ 0.32 which would correspond to a s [Mz) — 
0.118, close to the Standard value extracted from the Z peak physics. Here we adopt 
an alternative perspective. Varying a s (M T ) in a generous interval 0.25 to 0.36 we 
would find for the duality-violating component R T 

- 0.06 < -^SRt < 0.07 . (8) 

We thus see that the violation of local duality at the r mass scale turns out rather 
small, below the 10% level. 

This experimental evidence does not guarantee that violation of local duality is 
universally suppressed by such a factor for all types of processes; on the contrary, it 
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is expected to be sensitive to the details of the process M, [TÖL O. Another considera 



tion is even more important. As will be stated later, effects of local duality violation 
must oscillate as a function of energy scale with vanishing averages. Therefore, an 
accidental suppression at a fixed mass scale just near the actual M T ~ 1.77 GeV 
cannot be excluded.0 

In spite of these caveats, the hadronic r decay width provides a rather direct 
confirmation that local duality must be a reasonable approximation above 2 GeV in 
similar circumstances, in accord with theoretical estimates [TOj ] . 



3 Executive Summary on Local Duality in Heavy 
Quarks 

At first it might seem that duality between the quark and hadron description has 
little chance to hold for total decay widths. For the OPE yields |Ï3| 




V{H Q )^m% 1 + -5- (9) 



- i.e., that the width of a heavy flavor hadron is controlled by the fifth power of 
the mass of the heavy quark. Yet the major part of this high power comes from the 
available phase space (the fact most obvious in the semileptonic decays Hq — > iv X) 
which is determined by hadronic masses, in particular by Mh q rather than m,Q. 
Nevertheless it turns out that summation over different hadronic decay channels 
yields the widths computed at the quark level with their nonperturbative corrections 
order by order in 1/tuq (more accurately, in the inverse energy release). This is due 
to a conspiracy between the strong interaction effects in the decay amplitudes and 
the hadron masses, which can explicitly be traced in semileptonic decays |T3] where 
the conspiracy is enforced by heavy quark sum rules |Ï5|| . 

As explained in the Introduction, the crucial criterion for the theoretical analysis 
is whether a reaction can be treated by the OPE. The latter does apply to sufficiently 
inclusive decay widths - be they semileptonic or nonleptonic. In this context the 
OPE is meant to yield an expansion solely in terms of expectation values of local 
heavy quark operators evaluated for the actual hadron Hq. 

The practical implementation of the OPE expresses the widths as an expansion 
in 1/rriQ (or the inverse of the energy release E r ) with the coefficients shaped by 
short- dist ance physics accounted for perturbatively: 



oo / .. \k OO i 

k=0 m Q 1=0 7r 



3 Available cxperimental data allow to limit the magnitude of duality violation in R T (M) also 
at the masses Al somcwhat below than M T , although the limit suffers from more experimental 
uncertaintics. 
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where A is a generic (dimensionless) quantity and (fik) k are related to the nonper- 
turbative expectation values. 

With only a few terms known explicitly in these series, there are obvious limita- 
tions in the theoretical accuracy. This problem has little relevance to the peculiarities 
of local duality, and is not even specific to OPE-treatable observables as compared 
to a generic infrared-safe quantity Duality violations are effect s over and above 
this trivial reason. Even in an ideal scenario where all terms in both expansions 
of Eq. ( pro| ) were known, the OPE series would not define the exact result A(mg) 
completely; the power expansion of the "practical" OPE is only asymptotic, even in 
the Euclidean situation ||. ft is this ambiguous component that behaves differently 
for Euclidean and Minkowskian amplitudes. 

Let us list certain features of local duality that have been clarified by theoretical 
considerations over the last years: 

(i) The primary criterion for addressing duality violation is the existence of the 
OPE for the particular observable. 

(ii) In general local duality cannot be exact at finite masses. ft represents an 
approximation the accuracy of which will increase with the energy scales in a way 
that depends on the process in question. 

(iii) Effect s of violation of local duality can only have the form of an oscillating 
function of energy (uiq, E r , ...), or have to be exponentially suppressed. Duality 
violations cannot be blamed for a systematic excess or deficit in the decay rates. For 
example, no local duality violation can convert tuq into Mh q in the total width. 

(iv) The oscillating component violating local duality may be only power sup- 
pressed. In real QCD it nevertheless is to become exponentially suppressed as well 
at large enough energy, fading out as e~\ E / M V with a positive 7. The onset of that 
regime, however, can be larger than the typical hadronic scale ~ 1 GeV - for exam- 
ple, it may grow with increasing N c . The details of the asymptotics, in particular 
the power of energy by which duality violation is suppressed, depends on underlying 
strong dynamics and on the concrete process. This power is rather high in total 
semileptonic widths of heavy quarks. 

(v) The OPE equally applies to semileptonic and nonleptonic total decay rates. 
Likewise, both widths are subject to violation of local duality. The difference here 
is quantitative rather than qualitative; at finite heavy quark masses corrections are 
generally larger in the nonleptonic decays. In particular, local duality violation in 
nonleptonic decays can be boosted by the accidental near-by presence of a narrow 
hadronic resonance. Similar effects are extremely suppressed for semileptonic decays. 

(vi) It is not necessary to have a proliferation of decay channels to reach the onset 
of duality, either approximate or asymptotic. An instructive example is provided 



by the so-called Small Velocity (SV) kinematics in the semileptonic decays fÏ6| . A 
complementary nonleptonic example was identified in the exactly solvable 't Hooft 
model jÏ7|. 

(vii) A divergence in the power expansion of "practical" OPE underlying vi- 
olations of local duality is related to the presence of singularities in the quark or 
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gluon propagators at finite (or even infinitely large) space-time intervals.^ This is 
in contrast to finite-order OPE terms which account for the singularities of interac- 
tions (for the perturbative corrections), or for the expansion of propagators near zero 
space-time intervals. A certain class of nonperturbative effects (presumably strongly 
suppressed) comes from small-size instantons, which are neglected in the simplest 
version of the "practical" OPE. They are not specific to local duality and are similar 
in both Minkowskian and Euclidean amplitudes. Such instantons contribute to the 
Wilson coefficients in the OPE computed in the short-distance expansion. 

(viii) The presence of such fmite-distance singularities is a general, rather than 
exceptional feature of theories with nontrivial self-interaction. In particular, they 
are not directly related to quark-gluon confmement. Consequently, duality violation 
exists in general even without confmement. The latter, however, may essentially 
influence the nature of the singularities. 

With the number of constraints which have to be satisfied by quark-hadron 
duality, the magnitude of the effects in r s i(& — > c) attributable to local duality 
violation proper can hardly reach the half percent level, and is expected even much 
smaller in r s i(fe— >u). The uncertainties associated with the general quark-hadron 
duality turn out much larger for exclusive B — > amplitudes due to the charm 
quark being only marginally heavy. 

Some of these fmdings will be explained in the following sections. 

4 Theoretical tools 

The central tooi in describing heavy flavor decays is the OPE yielding series in pow- 
ers of 1/rriQ. The Standard application of the OPE is expanding deeply Euclidean 
Green functions at small space-time intervals where it yields well defined asymp- 
totic series. On the other hand decay amplitudes are shaped by processes in actual 
Minkowski space, and in general are determined by long-distance dynamics inde- 
pendent of the quark mass scale. Expansions in the heavy quark mass yield only 
very limited constraints under such circumstances. Yet inclusive heavy flavor tran- 
sitions are controlled by short distance processes, albeit in Minkowski space. The 
Minkowskian domain of physical interest and the Euclidean regime where the OPE 
is well formulated are connected through dispersion relations based on the analytic 
properties of the transition amplitudes governing the inclusive decay probabilities. 

To describe generic inclusive rates one starts out by defining a transition oper- 
ator T through the time ordered product of the appropriate weak Lagrangians L w 
responsible for the selected class of weak decays 

T{üj) = I d 4 xe-^T{L4x),Lt(0)} (11) 

4 The divergence has nothing to do with the ultraviolet divergence of separate higher-ordcr 
condensates evaluated without an ultraviolet cutoff, but is rootcd in their factorial growth with 
the order in the 1 /rriQ expansion. 
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and the corresponding forward "scattering amplitude" off the decaying hadron at 
rest 

A ^) = ^-{H Q \f{u)\H Q ) . (12) 
zm h q 

The auxiliary variable uo is defined in the complex plane and is the counterpart of 
the energy in r decays. In some respects it is similar to a variable mass of the heavy 
quark. A{u) has usual analytic properties in the complex u> plane. lts discontinuity 
at the physical point uj = gives the inclusive decay rate: 

r ** = 7^E( 2 ^ 6 ^ P n-PH Q )(H Q \LU0)\n)(n\L w (0)\H Q ) = 2W(0) , (13) 

ZM H Q n 

where n generically denotes all final states. 

On the other hand, the OPE expands A(ou) in powers of l/(u — E T ) with E r 
denoting the energy release in the weak decay. As mentioned above, such expansions 
are well behaved at complex ui away from physical cuts. Alternatively, A(uj) can be 
expressed as a dispersion integral over its cuts where the discontinuity is given by 
the total probabilities of the weak processes. This yields the power expansion for 
the inclusive decay rate. 

The local duality violations emerge as specific components of the decay widths as 
functions of energy, to which the dispersion integrals are poorly sensitive. In partic- 
ular, they can be due to the singularities in the forward "scattering" amplitudes near 
the physical point ui = 0|] that is near the masses trq where the decay amplitudes 
are singular, due to the thresholds for new channels or narrow resonances. 

In nonleptonic decays the operator T(u) is most generally given by the product 
of three Green functions of the final-state quarks, contracted with Q(x) and Q(0). 
The quark propagators are affected by a gluon medium or by quark condensates, 
Fig. 1. The singularities at physical uj can a priori be quite arbitrary and essentially 
depend on strong dynamics. For example, A{uj) can have a pole l/(rriQ— m + |T) if 
there is an appropriate resonance at the mass tjiq. Such a contribution is accounted 
for by the OPE which, in practice effectively averages the width. Therefore, the 
OPE does not distinguish a regular resonance structure from a smooth behavior 
obtained by averaging the width over an interval of rriQ. As a result, the deviation 
of the nonleptonic widths from the OPE predictions at certain discrete values of mg 
can, in principle, be arbitrarily large if narrow resonances existed with large enough 
masses. 

The situation is quite different in the semileptonic widthsF] where two of the 
Green functions are replaced by free lepton propagators, in the coordinate 

representation. This essentially limits the structure of the singularities and strongly 

5 The variable uj can be visualized as a shift in the decaying quark mass, toq — * tuq — uj; this is 
an approximate identification affected by 1/mg corrections. Therefore, position of singularities at 
real uj determincs heavy quark masses where nonanalytic behavior occurs. 

6 Or for decays like b — > s + 7. b — > s + U which we do not address here. 
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X X 

Figure 1: Examplcs of the diagrams generating the transition operator T(ui). The variable u> 
can be thought to enter in the weak vertices where the amount of (complex) cnergy is sucked away 
or pumped in. Computing the inclusive width we take the absorptive part of the corresponding 
amplitude. 



suppresses them at physical 10. The maximal possible violation of local duality is 
drastically reduced. 

4.1 Semileptonic decays — general formalism 

In semileptonic transitions the lepton momentum is an additional kinematic variable. 
In order to study the differential distributions, one considers a simpler forward 
scattering amplitude involving two weak quark currents: 

V = ^(#d j^xe-^iT{j,{x)Jl{0)}\H Q ), (14) 

with = ^7^(1— 7ö)<5 f° r Q ~^ Q weak decays. From now on we assume explicitly 
the case of B meson decays, and also neglect, for simplicity, lepton masses. Decays 
leading to a tv t pair do not introducé essential new features. 

The tensor h^ v is generally decomposed into five invariant functions hi depending 
on k 2 and ko = kv: 

hfu, = -hig^ + h 2 v il v v - ih 3 e lu/a pv a k 13 + h A k^k v + h 5 (k^v u + v^k u ) , (15) 

with denoting the 4-velocity of the decaying meson, and their discontinuities are 
the heavy quark structure functions: 

Wi — 2Im/tj . (16) 

For massless leptons u>4 and w§ do not enter; moreover, W3 does not affect the decay 
width: 



G 2 F | V q b | 2 



r sl = *>":' dk 



.2 



dk yjk$-k 2 (k 2 Wl (k ] k 2 ) + \(k 2 -k 2 )w 2 (k ; k 2 ) 



167T 4 Jk 2 >0 Jk >Vk? V 3 

n? WA 2 r°° _ _ r / -, -, 1 

(J F\ V qb\ / aZ1\Z\ l AU I/U2 £2\.. / 7 „ . 7„2 t2\ , L 



[°°dk 2 \k \ [ _ dkJ(k 2 -k 2 )w 1 (k ;k 2 -k 2 ) + -k 2 w 2 (k ;k 2 -k 2 )) 

JO Jk >\k\ \ 6 J 

(17) 

(the upper limits of integrations are fixed by dynamics rather than by kinematic 



167T 4 JO Jk Q >\k\ 

ons are fixed bj 

constraints; namely, Wi vanish for k > (M| + Mp — k 2 )/2M B , and likewise for 
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b — ► u.) This expression shows that the singularities in the semileptonic width are 
suppressed compared to what one finds in the general case for nonleptonic decays. 
The contribution of the particular hadronic state X with the invariant mass Mx to 
the structure functions w(ko] k 2 ) has the form 

2 \ I 2 r / i M%-M x + e\ 



5W(k ;k 2 ) oc \F x (k z )\ 2 5 k - 



2M F 



where Fx denote the B — > X transition amplitudes up to certain kinematic factors. 
Since Fx(k 2 ) are analytic at physical k 2 accessible in the decays, the integration in 
Eqs. ([Ï7| ) provides an amount of smearing suppressing the threshold singularities. 

Representation (|ïj]) itself is not, however, most convenient to reveal the related 
suppression, since the interval of integration over ko shrinks to zero as k 2 approaches 
E 2 = {m h — m c ) 2 . The duality properties of the total semileptonic width are more 
explicit if it is represented in the following form: 

r sl = ƒ dk dk 2 ^k 2 (k 2 Wl (k ; k 2 ) + Ukl-k 2 )w 2 {k Q] e 

rfl \y ,|2 F kl _ / i _ 

U F\ v qb\ / AU / A t2\u\ I fl.2 ÏT2\„„ /r. . i.2 £2\ , 1 £2 „„ /r, . r.2 7> 



16tH 



/ d£; / dk 2 \k\({k 2 -k 2 )w 1 (k ;k 2 -k 2 ) + -k 2 w 2 {k ;k 2 -k'' 
Jo jo V 3 



(19) 



or 



G \V h\ f°° 

r sl = Ta 4 / dfc U ( k o) , 

ion Jo 

u(k ) = dk 2 ^k 2 -k 2 (k 2 Wl (k 0] k 2 ) + l(k 2 -k 2 )w 2 (k 0] k 2 ) 

dP|fc|((fc 2 -P) Wl (fc ;^P) + ^ 2 «; 2 (fco;^-fc 2 )) . (20) 



As in the Standard representation (|Ï7|) , the hadronic function u(ko) actually vanishes 
at ko > Mb — Md. Likewise it can be represented as a discontinuity of a function 
W(ko) related to the transition amplitude: 

W(k ) = £° dk 2 ^k 2 -k 2 (k'hiko; k 2 ) + ^(k 2 -k 2 )h 2 (k ; k 2 )^j , (21) 

u(k ) = - \w.(k )-W + (k ) 

where W±(ko) are obtained by integrating over ko above and below the real axis, 
respectively. W(k ) has analytic properties similar to the transition amplitude h^ u 
at fixed k 2 or k 2 . In particular, since the integration in Eq. (|2Ï|) has fixed end 
points, W(k ) acquires singularities at k where h(k , 0) or h(k , k^) become singular. 
The former corresponds to the thresholds with vanishing lepton invariant mass k 2 : 
k = x ■ The l^ter corresponds to the zero-recoil thresholds ko = Mb — Mx- 
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It actually has a special explicit counterpart in the OPE since it represents the 
maximal-/c 2 part of the decay probability [D|. However, due to the integration over 
k in the first of Eqs. (pÜ|), this singularity is not important, and we do not go into 
further details here. 

The total width is obtained as an integral over the discontinuity u(ko) of W(ko) 
along the cut down to k = 0, and therefore can be represented as a contour integral 



G 2 p\ V q b\ A 1 



16tH 



- / dk W(k 



[22] 



c 



see Fig. 2. This integral is not an analytic function of energy release or of the heavy 
quark masses since the final points of the contour are fixed at k = 0±z0 lying on the 
cut. However, this shows that the singularities of the width and, correspondingly 
the violation of local duality is determined by the transitions with k 2 = (say, the 
thresholds with Mx = Mb) in spite of the kinematics with k 2 > contributing to the 
total rate. (A similar fact was explicitly observed for instanton effects in Ref. ||.) 
The full hadronic energy release m^ — m c controls the total width. 



C 



Figure 2: The integration contour around the cut of W{ko) for the total semileptonic width. 
The end points at ka = are shiftcd into the complex plane for clarity. 



As before, this can also be shown by analytically continuing the correlator asso- 
ciated with the total width, using the auxiliary variable to. The amplitude A{oj) of 
Eq. ( |Ï"2"D is literally given by 



(23) 



A(lu) = -9ËYé^ [ dk dk 2 \k\\n(k 2 -k 2 ) Ukl-k^h^+uj-kl-k 2 ) + 

^k 2 h 2 (k +uj; k\-k 2 \ 
For our purposes it can instead be defined directly as 

A ^ = -^TJ d H dk2 "J k °- k2 {k 2 h(k +u;; k 2 ) + -(k 2 -k 2 )h 2 (k +u;; k 2 )J 



or 



Cl 2 \V A 2 r M r( k o^) 2 i / 

A(u) = / dA; / dk 2 J(k -uj) 2 -k 2 (k 2 h(k ; k 2 ) + 



^(k -u) 2 -k 2 )h 2 (k ;k 2 )^ (24) 
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(M is a sufficiently large constant), so that for real u 



2 \mA{u) = -^-f- / dk ü(k ; W ) (25) 

with 

«(*<,;«;) = ^ dfc 2 v /(fc -ü;) 2 -A; 2 ( fcV^ol fc 2 ) + -((A; -^) 2 - A; 2 )^ 2 (A; ; A; 2 )j . 

(26) 

The approximate scaling holds «4.(a>;m&) — A(0;mb — uj;), and one obviously has 
ü(k ; 0) = u(k ). 

4.2 Sum rules 

Sum rules for semileptonic transition amplitudes refer to a somewhat different type 
of inclusive probabilities sharing some peculiarities of total semileptonic widths. 
Besides providing a wealth of important dynamic constraints, they help to illustrate 
and visualize a number of the OPE results obtained directly for the total widths. 
In particular, they allowed to show in the most general setting that the total decay 
probability coincides with that of the free quark ansatz, and how the hadron masses 
get replaced by the quark masses thus extending this equality to the level of I/tuq 
corrections. 

In sum rules we fix one of the two kinematic variables, usually k 2 or k 2 , and then 
consider the total transition probability weighted with a power of energy ko, say 

Wik 2 ) = [dk (-k o r Wl (k ; kl-k 2 ) (27) 

Z7T J 

(in the heavy quark expansion we actually count energy from a different value rather 

than ko = 0: we put e = m^—^m 2 ^ + k 2 —ko so that e = corresponds to the free quark 
kinematics). Dispersion relations directly equate these moments to the coefficients 
of the asymptotic expansion of the hadronic functions hi in powers of 1/fco at large 
(complex) k (or e) computed in the OPE. The moments considered in the sum 
rules are similar to inclusive differential decay widths at fixed k 2 (or k 2 ). However, 
they are not restricted by the kinematic constraints and include states with high 
mass which may not be accessible in the actual decay, but only in the scattering 
of the weak current on the heavy hadron. For this reason, the sum rules are in a 
sense exact relations not affected by local duality violations. Yet physics of duality 
is present in the sum rules as a question at which mass scale and how accurately 
the integrals are saturated and where their normalization point dependence enters 
the perturbative stage. This scale is directly related to the onset of local duality in 
the actual total semileptonic widths. 

Another aspect of duality in the sum rules, again not directly related to local du- 
ality, is of how accurate is the expansion in the inverse charm mass. We can compute 
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the moments in QCD expanding the transition operator in powers of KqcY>/E q with 



E g = ym q + k 2 ; in practice one typically has E q ~ m c . Here one encounters a new 
aspect called global duality [15]. It was shown that the quark-hadron correspondence 
holds separately for the decay-type charm intermediate states and those which are 
due to other physical processes in the full transition amplitude h^ u , order by order 
in the 1/E q expansion. However, the identification may be lost in the exponential 
terms cx exp (-2£ 9 //i hadr ). 



4.3 Small Velocity expansion in b — > c transitions 



Many simplifications arise for b — > c semileptonic decays when the velocity of the 
final state hadronic system is small with both initial and final state quarks sufficiently 
heavy. The physics of these SV transitions is also more transparent | Ï6| . The SV 
case is a particularly relevant testground for duality violation: the OPE for total 
inclusive widths relies on the expansion in the energy release m&— m c rather than in 
m,b itself. In addition, the peculiarities of the four-fermion decay interaction which 
manifest themselves in the high power n = 5 of the dependence on the quark masses, 
suggest ]Ï4| that in general the actual hardness of the total widths can be further 
decreased down to only a fraction of mb- This a priori can boost effect s of local 
duality violation, but by the same token enforces the onset of the SV kinematics. 
Here we briefly describe the SV regime for the total semileptonic width. 

The heavy quark expansion of the semileptonic b —* c width is usually written 
as an expansion in l/m&: 



r s i(-B — » x c 
T (b - c) 



iA 



c (m c /m h ) + c 2 (m c /m b ) — ^ + c 3 (m c /m b ) 



mi 



i4 

ml 



+ 



(28) 



where the coëfficiënt functions c& are computed in the short-distance expansion 
depending on a s , and on m c /mb as a parameter. There is no contribution of order 
1/mb - the principal result of the OPE [0. In the SV expansion we rearrange terms 
in Eq. (28) using instead the parameters E r = mb — m c and rrib', this implies that 
E r rrib holds in the applications. The velocity of the final state hadrons does not 
exceed (mb — m c )/m c ~ E r /rrib 1, and we have 



,(!) 



El 



v. 



(2) 



E r 



El 



m b -m c 



V. 



(3) „(3) 



+ 



El El 



(29) 



m b m b 

Here ujp are hadronic parameters of dimension Aq CD . The OPE relates them to 
the expectation values of local heavy quark operators in the B meson state. This 
expansion is the nonrelativistic version of the OPE of Refs. [ |T5| , [TJ5j kinematically 
similar to the neutron /3-decay, improved by accounting for the relativistic effects 



16 



order by order in velocity. The SV semileptonic width was considered to order v 2 by 
Le Yaouanc et al. [2Ö| in a toy model where the OPE generalities become explicit. 

Expansion (p9|) exhibits a few peculiar features. First, there are no nonpertur- 
bative corrections in the velocity independent term - i.e., even at mj- m c — > 0, 
for large enough heavy quark mass. The impact of strong dynamics here is given 
solely by the short- dist ance factor N = 1 + \g\ with g& denoting the perturbative 
renormalization of the (diagonal) axial current. 

Secondly, there is a unique term linear in v; it is related to the quark spin and 
is absent for spinless quarks (or for vector-like weak interactions) . 

The OPE explicitly performed for T s \ through order 1/mf gives concrete expres- 
sions for the terms ~ 1/-Ey and ~ 1/Ef, with arbitrary powers of v via the kinetic, 
chromomagnetic and Darwin expectation values in B meson. The SV OPE also 
exhibits distinct features. Namely, the OPE series to a given order in v have only a 
finite number of terms[| in 1/E r , although this number grows with the power of v. 
In particular, there is a single term to order v 1 as anticipated in Eq. (|29|), which is 
given by /i 2 G : 

(1) 10# A 2 / 2\ fon\ 

"2 = ~ i +3g 2 M m 6 ) (30) 

(there are already seven terms through order v 2 ). 

Returning to violations of local duality in the semileptonic decay widths, it is 
likewise advantageous to analyze it as the difference between the exact width and 
its OPE expansion, both considered to a given order in v. This has an added 
convenience since then one does not have to address the general problem of defining 
an exact sum of the asymptotic expansion with factorially growing coefficients 
In this respect, for any particular term ~ v l the situation is similar to the case of the 
r decay width considered earlier. Similarly, the fact is evident that duality violation 
is to be present at some level regardless of details of dynamics.0 

As mentioned above, duality is exact to order v°. It is less trivial a priori 
that the same holds even to the next order in v: There is no difference between 
r s i and Tope up to terms v 2 at any energy release! This can be seen by noting 
that only the "quasielastic" widths B — ► D and B — > D* must be considered. 
The chromomagnetic interaction shifts masses in the same direction for B and D; 
therefore, this effect shows up only to ö(v 2 ). However, S^Mp* = —^^S^Mb, and 
this splitting is the only source of ö(v) corrections. (Recall that the phase space is 
controlled by the fifth power of E r = Exactly the same contribution evidently 

emerges in the OPE to order v 1 . 

The actual violation of local duality in F B \ emerges only to order v 2 and is 
suppressed by at least two powers of E r . Since = ^ these effect s start with 

7 This follows from the second representation in Eq. jÏ9|). In the SV expansion the structure 
functions are polynomials in k; integrating them over k yields a finite polynomial in fco whosc 
nonperturbative expansion therefore has only a limited number of terms. 

8 The asymptotic nature of the expansion in velocity Eq. (|2^) itself is an additional source of 
duality violation. Since it is not instructive, we do not discuss it here. 
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terms oc 1/mf; higher order terms can be further suppressed already by factors 
A i hadr/( m b — m c)- in actual B decays E r ~ 3.5 GeV, i.e. much higher than, say, the r 
mass; this underlies the strong numerical suppression of the local duality violation 
in Y A {B). 

To order v 2 the total decay rate is given by a few moments of the nonrelativistic 
SV structure functions. (The exact combination of moments is determined by the 
nonrelativistic form of Eqs. (0).) They are computed, however, over a limited 
energy range. If the energy release exceeds the mass of the heaviest charm state 
(the "P wave" excitations contribute here), these moments exactly coincide with 
their OPE expressions. In practice, there are always high enough charm states, for 
instance, dual to perturbative gluon excitations of charm. Therefore, the onset of 
local duality in the SV width is directly related to the scale where nonperturbative 
contributions to the sum rules are saturated and to onset of the perturbative regime 
there, the dashed-dotted line in Fig. 3. We will refer to this figure more than once 
later, and its meaning will be explained in more detail there. 

5 Theoretical implementations 

As illustrated above, while the concept of quark-hadron duality can be formulated 
accurately whenever the OPE can be applied, it represents a complex phenomenon. 
Studying it in a concrete model would be of great help. However this represents 
a nontrivial task; for some aspects of the OPE and thus also of duality depend on 
subtle features rooted in the gauge nature of QCD. Since those are typically not 
respected in naive quark models, doubts arise in the relevance of their numerical es- 
timates of duality violations. Below we describe two other dynamic implementations 
of nonperturbative physics. While both are highly nontrivial, neither contains the 
full complexity of real QCD. Yet they are largely complementary in accounting for 
different facets of nonperturbative QCD. More details can be found in the already 
mentioned review by Shifman M or in the original publications referred to there. 

5.1 Instanton model 

Instantons provide a nontrivial dynamic realization of nonperturbative physics gen- 
erating full OPE series which, in principle, can be evaluated to sufficiently high 
order. They are believed to be relevant for the properties of low and intermediate 
energy physics shaping a host of nonperturbative parameters of QCD in Euclidean 
space. Models of this type assume that quarks and gluons propagate not in the 
perturbative vacuüm, but in the background of instanton configurations of typical 
size po- The instanton density decreases fast for instantons of a smaller size. It is 
usually assumed that the density of instantons is small enough compared to 1/pg 
in order to preserve a meaningful notion of individual instantons or to facilitate the 
computations. In this approximation the instanton effects are proportional to their 
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density and strongly depend on their mean radius po- 

The instanton ansatz is instructive in a number of aspects related to duality and 
its violations [§]. 

(i) It has been explicitly shown to lead to duality violations in total semileptonic 
widths contrary to the often stated lore that while nonleptonic widths suffer from 
duality violations, semileptonic ones do not. The instanton calculus has demon- 
strated a quantitative rather than qualitative distinction in this respect between the 
two processes. 

(ii) The instanton ansatz illustrated that fmite-distance singularities in the Green 
functions - leading to divergences of the OPE and to violations of local duality - are a 
common, rather than exceptional feature of strongly interacting systems. Moreover, 
it showed that duality violations are not intrinsically tied to confinement, contrary 
to what a historical perspective might suggest. At the same time, the instanton 
ansatz exhibits the general features discussed above - oscillations and fast decrease 
with energy, strong suppression upon averaging, larger effects in nonleptonic widths 
compared to the semileptonic ones, etc. 

(iii) On the practical side, it was shown that conventional instantons cannot 
induce any appreciable duality violating effects in total semileptonic widths of B 
particles, regardless of uncertainties in the model parameters. Even boosting up the 
possible instanton density leaves the effects below the permill level. If some appre- 
ciable oscillation observed in e + e~ or r decay distributions are rooted in such effects, 
the responsible nonperturbative configurations must have significantly smaller size. 
While having a minor effect on Euclidean short- dist ance observables at energy scales 
around rrib, they would be manifest at a lower scale ~ m c even in the Euclidean do- 
main. 



5.2 't Hooft model 



The 't Hooft model [fJTJ - a (l + l)-dimensional analogue of QCD with N c — > oo - is 



an attractive theorist's laboratory in exploring various complicated aspects of non- 
perturbative dynamics in QCD. Being solvable, it allows in principle deriving pre- 
cise numerical values of the modei's counterparts of actual hadronic characteristics. 
Then it is possible to confront them with the results of particular approximations 
employed in real QCD and in this way to test their viability. 

This model is particularly appealing for studying duality and its limitations. It 
automatically respects the basic underlying features of QCD related to gauge invari- 
ance, including its rigorous sum rules; ad hoe models typically fail in that respect. 
At the same time, there is little 'wiggle room' for adjusting various parameters, so 
that the results are subject to smaller interpretational bias. 

Since the final states in the 't Hooft model consist of an infinite series of dis- 
crete narrow resonances, one expects this model actually to maximize local duality 
violations. Additional considerations can be found in the dedicated papers |TÖ|, [17) . 

Since it is solvable, one can determine inclusive widths by calculating the rates 
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for all available exclusive channels and sum over them. Comparing the sum with 
the OPE result provides a direct test of duality. 

This program has been performed for both nonleptonic and semileptonic decays 
in the above papers, and full consistency with the OPE has been found for both 
semileptonic and nonleptonic decays. In particular, the 1/toq terms forbidden by 
the OPE are absent from the total widths.0 

At the same time, violation of local duality allowed by the OPE is explicitly 
observed in the 't Hooft model. The difference between the actual width and its 
OPE expansion is oscillating and fades out fast with the increase in energy. The 
decrease happens even without including finite widths for the resonances or smearing 
in energy. The degree of suppression, however, is non-universal and depends on the 
type of the transition. 

The numerical aspects of duality in the semileptonic decays in the domain of 
intermediate masses were analyzed in Ref. p2| . It has been established that the 
prominent resonance dominance does not spoil duality, nor delays its onset. The 
key fact is the saturation of a few lowest sum rules: above this scale the amount 
of duality violation is negligible. Violation of local duality in the total semileptonic 
B decay widths was found to be safely below the permill level even considering the 
possible uncertainties in translating quark masses. It should be noted that in general 
the nonperturbative effects in the 't Hooft model are quite significant, in particular, 
in exclusive heavy quark amplitudes. Inclusive decay widths appear to belong to the 
class of observables most robust against the effects of nonperturbative dynamics. 



6 Application to T S \(B) 

A central phenomenological application of the heavy quark expansion in beauty 
decays is extracting the value of the CKM parameter V c t,. The total semileptonic 
width of B mesons is a well measured quantity that can be treated theoretically 
through the OPE. Thus it is important to scrutinize the possible size of local duality 
violation in F S \(B). 



6.1 General features 

While the exact mechanism driving violation of local duality in QCD is not reliably 
known, it has many features which significantly limit its magnitude in semileptonic 
decays. We illustrate this first semiquantitatively using the general arguments of 
Ref. []Ï2| , Sect. 2.5.3; in the next subsection we present a more specific estimate. 



The underlying expansion parameter for T S \(B) is the energy release E r = m^ — 
m c ~ 3.5 GeV rather than the safely large mj. Therefore, we can analyze the amount 

9 It had been claimed in the literature (and most recently reiterated by Grinstein) that the OPE 
is inapplicable to inclusive decay widths, based on the alleged presence of 1/mg ternis in the total 
widths within the 't Hooft model as deduced from some numerical computations. Thosc findings 
were erroneous. 
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of duality violation varying m c from to m;, while keeping Mg fixed to its actual 
value: 

XC^r) = ^ 77 7TT - 1 • (31) 

For m c — > 0, i.e. effectively for 6 — > w, the magnitude of duality violation is 
expected to be quite negligible, since the energy release is plentiful. We shall con- 
servatively assume it to be below 1CT 3 , even though all realistic estimates yield much 
smaller values || |2Ï| . The opposite side of the interval - m c — > - deserves care- 
ful attention, presumably yielding the largest values of \- One still has to maintain 
nih — nic > /ihadr to keep the literal OPE series meaningful.0 

For mb — m c ~ 1 GeV one has v ~ 0.2 for the velocity of the c quark and the 
SV expansion Eq. (^) becomes a valuable tooi. As explained in Sect. 4.3 duality is 
exact to orders ö(v°) and ö(v 1 ). It is worth noting that the leading nonperturbative 
term in Eq. (|29|) amounts to about 15% at m?, — m c = 1 GeV (5% at the actual 
mfe— m c = 3.5 GeV), and duality is still exact at this level! 

Violation of local duality as measured by \ emerges first at ö(v 2 ). The typical 
scale of the relevant nonperturbative effects can be estimated by taking just the 
leading OPE term. With 

4 2) - ~/4 - ~/4 02) 

it is around 2%. In reality, the duality violating component must be even further 
suppressed. 

Through order v 2 both the OPE terms and the individual decay probabilities are 
determined by the SV amplitudes - the transition formfactors called T3 and ri to 
the excited states in the heavy quark limit (the exact definitions can be found, for 
example, in the review ||23|| ). They have been well studied theoretically. The OPE 
parameters v\ , v\ are given by the corresponding moments of the SV structure 
functions, which are sums of terms like Ir^l 2 ^. The very same sums appear in the 
total width computed as a sum of the quasielastic [B — > D,D*) and P-wave decay 
probabilities. The only difference is that the actual decay width includes in the sum 
only the kinematically allowed states, those for which ^ m& — m c . Thus, if the 
SV sum rules were completely saturated at a given energy release, duality violation 
would have been totally absent at this level, and x ^ 0(v 3 ). 

The first P-wave states have excitation energy e\ ~ 400 to 500 MeV. In the 
't Hooft model, they almost completely saturate the sum rules which implies minute 
duality violation even at the minimal energy release. This probably does not hold 
for real QCD, where a significant contribution can be attributed to the states with 
e up to 700 to 800 MeV, in particular for the higher moments. Let us note that in 
traditional analyses (say, in the QCD sum rules) the sum rules are assumed to be 

A 2 

10 It is curious to note that even at Aqcd > TUf, — m c S> ^ p ~ 40 MeV duality holds up to 
terms ~ ; assuming fJ-f ladl ~ 0.5 GeV 2 one would estimate this effect as only a few percent. 
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saturated within ten percent by the excitations below 1 GeV. This would lead to a 
similar minute duality violation at the sub-10 -3 level already at m^ — m c ~ 1 GeV. 

We do not go that far; instead we assume v\ to be merely 75% saturated, and 
allow an even more modest saturation of 50% for higher moments. Then we arrivé 
at 

|x(lGeV)|£l% (33) 

A more refined bound can be obtained applying reasoning similar to the one we em- 
ploy below, Sect. 6.2, making use of general properties of duality violation: one can 
relate the maximal magnitude of the duality violation to a fraction of the branching 
fraction of the highest open threshold. Namely, to order v 2 we have 

X (E r ) = -^Q (2) (£ r ), (34) 

where - analogous to r decays - is an oscillating function with a known threshold 
behavior at each P wave state, and its first few moments vanish. Although this 
results in a more stringent bound, we do not present it here. 

The terms of higher order in the velocity v are suppressed from the very begin- 
ning, v 3 ~ 0.01 at mfc— m c = 1 GeV. Therefore, we conclude that the magnitude of 
duality violation is at most a fraction of a percent already at mb — m c = 1 GeV. 

Since the amplitude of duality violating oscillations decreases fast with energy 
release, we see why it always emerges extremely small for b — > u transitions. With 
the amplitudes decreasing as a power of energy release, the scale of duality violation 
at physical masses, E r ~ 3.5 GeV is found well below a phenomenologically relevant 
level. Some additional qualitative arguments including a comparison with r decay 



width can be found in Ref. [12 



6.1.1 Nonperturbative hadronic scale in semileptonic B widths and heavy 
quark symmetry in charm 

Fig. 3 will help to make our reasoning more transparent: it shows salient features 
of the typical hadronic mass distribution in B — > X c Iv. 

The canonical assumption is that the energy scales above 1 GeV beyond Md» 
belong to the perturbative regime. Yet based on the QCD Lagrangian alone we 
could not rule out a priori that, for example, a prominent charm resonance existed 
with a mass above 4 GeV that is not shadowed by a depletion in the continuüm 
hadronic mass spectrum nearby, and so altogether the total yield were not dual 
to the perturbative spectrum. If such an unorthodox resonance exceeded Mb - or 
were at least close to it - this could lead to a poor convergence of the OPE and/or 
sizeable duality violation in T S \(B). In the language of Fig. 3 it would mean that 
the resonance region R extended up to or even beyond Mb with little or no room 
for the perturbative domain P. 

Based on our whole experience with QCD this is not a likely scenario, but its 
consequences can be analyzed. Such a phenomenon would leave heavy footprints in 
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other transitions involving charm as well. In particular it would drastically impair 
the calculability of the B — > D* formfactor even at zero recoil. Oddly enough, this 
connection is mostly missed in the literature. 

The point is that even the Euclidean heavy quark expansions depend on the 
same nonperturbative operators, and their expectation values are correlated with 
the spectrum and the residues of the heavy flavor resonances. Suppose we have 
a charm state X with Mx — 4.5 GeV and the amplitude (X|c7 M (l — ^y 5 )b\B) such 
that its branching fraction constitutes, say, only a minute 0.5% part of r s i(P). The 
velocity of such a hadron being very small, it must be either an 5"-wave or P-wave 
transition. In the heavy quark limit only P-wave amplitudes vanishing at threshold 
survive, 

J^(X(v)\c lfl (l- l5 )b(0)\B) kt x v. (35) 

However, there are l/m c corrections which yield the amplitude nonzero even at van- 
ishing recoil, and they are of order unity at such M x . Using the nonrelativistic 



expansion of the zero- recoil cb current one finds (see [15], Sect. VI) that such ampli- 
tudes are generally given in terms of Mx ^ mc r, up to spin-related factors. Therefore, 
we can simply assume the amplitude a constant at small v, and approximately equate 
it to the corresponding SV amplitude r in the heavy quark limit. The partial decay 
width can then be estimated as 

T sl (P - X) « \r\\M B -M x f , (36) 

that is, say, even for Mb~Mx — 1 GeV we would need 

r 2 x * 0.2 . (37) 

However, let us look now at the contribution of such nonperturbative states to 
the sum rules, see, e.g. [23[. The increase in the IW slope through the Bjorken sum 



rule is given by 

Ö X Q 2 « 3t x « 0.6 (38) 

and is, in principle, tolerable keeping in mind the tentative nature of the estimates. 
However, everything is different already for A = Mb—vti^ and, in particular, for /x 2 : 



Ó X A ~6r x e x « 3 GeV 
öxfil ~9r x e x «11 GeV 2 



(39) 



not to mention the Darwin term and other higher-order operators. On the other 
hand, the literal value of P(0) in the 1/m expansion shifts at least by — 1% for an 
increase in /z 2 by 0.12 GeV 2 (this is a model-independent bound f24|). It is then 



obvious that through the exact dispersion relations and the OPE the expansion of 
the B — > D* amplitude is orders of magnitude more sensitive to the high-momentum 
nonperturbative dynamics than the potential duality violation in the total semilep- 
tonic width. 
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The above estimates can be cast into a more accurate form using the whole set 



of the heavy quark sum rules including new exact sum rules [25|. We have partially 



incorporated the latter in Eqs. (|3q), (p9|) and (p6|) assigning approximately equal 
values to the corresponding ri and T3 for the highly excited states. If this were 

2 2 

violated, one would completely destroy the exact relation for the B meson spin and 
got an unacceptedly large value for the chromomagnetic operator \j? G . However, in 
view of the obvious trend of numbers, such an improvement seems superfluous. 

The key point illustrated by the above consideration is transparent. Physics of 
power corrections in individual b — > c semileptonic transitions is strongly correlated 
with the saturation of the heavy quark sum rules, since the former include the 
expectation values given, for example, by the moments of the SV structure functions. 
Allowing appreciable genuine nonperturbative effects in the heavy quark hadronic 
system with excitation energies /i exceeding a couple GeV would dramatically upset 
the l/m c expansion leading to the higher order terms scaling as powers of fj,/m c 
instead of the naive Aqcd/^o while only moderately affecting the total semileptonic 
B width. This simple interrelation was put forward already in the review f26j, but 
seems to be missed up to now by most of the heavy quark theorists speculating 
about duality violation in the total semileptonic widths. 



6.2 A realistic estimate of duality violations in T S \(B) 

The above considerations were intended to illustrate the scale of rather model- 
independent factors suppressing violations of local duality in inclusive semileptonic 
widths. Augmenting them with explicit estimates of the power terms and the asymp- 
totics of the nonperturbative formfactors yields an actually much stronger suppres- 
sion. The main lesson here parallels what has been derived in the 't Hooft model: 
the usual resonance-related duality violation becomes negligible as soon as the char- 
acteristic scale of their mass gap is passed JEJ. This scale /^hadr is believed to be 
about 0.7tolGeV in QCD. 

Let us take a look at Fig. 3: p(Mx c ) corresponds to the structure functions 
w of the b quark inside the B meson appropriately integrated over the spacelike 
momentum. It is related to the total decay width as a function of the energy 
release with the lepton phase space factored out. For clarity, we actually multiplied 
this quantity by the factor Mx c — rn c to make it approximately constant in the 
perturbative regime. The domain R below Md* + Aihadr (i-e., E r ^ /ihadr) is usually 
referred to as the resonance region, while the one above it is viewed as described by 
perturbative dynamics. 

Our preceding discussion concentrated on the duality violation which is asso- 
ciated with the resonance domain R. In a sense, we implied that the perturba- 
tive contributions were absent, and the distribution in Fig. 3 nearly vanishes above 
Md* + Athadr; thi s is similar to the situation in the 't Hooft model where the pro- 
duction of highly excited states with masses above m c + Aqcd is strongly power 
suppressed, and the principal nonperturbative effects originate from and can be re- 
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Figure 3: Plot illustrating different strong interaction domains. p(Mx a ) stands for a generic b — > 
c structure function of a beauty hadron, related to the observable distribution over hadronic final 
state mass up to kinematic factors. R marks the 'standard' nonperturbative domain dominated 
by low-lying resonances. In the quark models the distribution fades out very fast above it in the P 
domain (dash-dotted tail), but is larger in actual QCD where it is dual to the parton processes with 
hard gluons. The decay distribution there may also oscillate around the perturbative prediction 
shown by the dash line, up to higher energy scale, with a larger interval of local duality. The 
cxcitation energy Mx c —M D (,) at the borderline between the two domains is generically denotcd 

Mhadr- 



lated to a few lowest resonance states in the domain R. If the resonances were 
totally absent above the dash-dotted borderline Md* + /ihadr, the actual hadronic 
width would exactly coincide with its OPE expansion for E r > /i na dr- The numerical 
studies of the duality in the 't Hooft model p2| were undertaken just to elucidate 
how the prominent resonance structure in the nonperturbative domain would affect 
the convergence of the 'practical' OPE as we approach this borderline and if it could 
delay the onset of accurate local duality. 

In real QCD, however, the relevant excitation probabilities decrease not powerlike 
as in simple models, but much slower due to the perturbative effect s associated with 
the emission or exchange of hard gluons, as shown explicitly in Fig. 3. Physics of 
such processes has not been included so far. In actual B decays we are safely in 
the perturbative domain, as in Fig. 3, if we place Mb close the right edge of the 



plot. As suggested in Refs. JT^, [22] ) an appreciable violation of local duality could 
possibly emerge from this type of hadronic excitations. We present its estimate in 
this section. 

The instanton effect s addressed in Ref. would likewise populate the domain 
above the resonance region. However, they lead to oscillations which decrease very 
fast in magnitude and can be neglected. They would not be visible in the relevant 
part of the distribution in Fig. 3. 

The perturbative corrections per se do not vitiate the OPE even though they 
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generally lead to divergence of the spectral functions moments. The short- dist ance 
expansion, instead of being regular order by order at large energies, acquires loga- 
rithmic factors or, more generally, non-integer powers of energy yielding continuous 
imaginary parts. It does not generate any duality violation by itself as is most ob- 
viously illustrated by finite-order perturbation theory: with the parton processes 
yielding smooth continuous spectra of states, there is always exact duality at any 
energy, regardless of the presence of kinematically allowed and forbidden states. 
Since all the decay probabilities are smooth here, they coincide with their smeared 
averages, and there is no room for violation of local duality. The spectral density 
p(Mx c ) as it comes out perturbatively is shown in Fig. 3 by the dashed line. 

However, the excitations of gluonic degrees of freedom may, in principle, exhibit 
a certain resonance structure. If the typical energy scale for unfreezing gluonic 
degrees of freedom is well below the energy release, our general discussion in the 
previous sections applies, and duality violation can be neglected. On the other 
hand, the mass gap for the transverse gluon degrees of freedom can be larger than 
in the qq channels. Then, as conjectured in Ref. |Ï7), one would observe two scales 
in the onset of local duality: it first sets in rather early, however only approximately 
with a typical accuracy ~ «s/V ~ 10%; one needs to ascend to higher energies 
characterized by the "gluonic" mass gap to achieve full duality. 

In this scenario the resonance structure can be observed at a suppressed rate 
up to higher energies, and its average is dual to the probabilities of the gluon 
bremsstrahlung, see Fig. 3. In semileptonic heavy quark decays the latter rates 
decrease slowly with energy compared to a powerlike suppression of exciting the 
valence quark- antiquark mesons. Therefore, this would dominate the local dual- 
ity violation at intermediate energies. Here we present an estimate of the possible 
magnitude of these effects; the resulting upper bound is to a large extent model- 
independent. To maximize the possible effect, we replace the continuous oscillating 
behavior by a comb of 5-functions representing infinitely narrow resonance; however, 
upon average they reproduce the perturbative distribution. 

The basic idea for this estimate goes back to Ref. |TÜ[ where it was applied to 
the t decay width. It makes use of the following facts: 

• The average of Tb — Tope as a function of mass vanishes. 

• Tope is a smooth function given by a series in I/tuq, while Tb has known threshold 
singularities. Each partial decay width is a smooth function above threshold. 

• The principal thresholds are assumed equally spaced at large masses, and the 
corresponding transition amplitudes follow a power-like scaling. Then rfl r T^ i ° PE is, 
to leading order in 1/wiq, a periodic function up to an overall suppression factor 
1/rriQ with some index k > 0. 

This allows to determine the actual leading order behavior of the duality violation. 
The index k is given by the (maximal) rate of the last open channel relative to the 
total decay rate - it also scales like I/tuq. 

It should be noted that the exact expression for Tope generally requires addi- 
tional definition since the series are asymptotic. The best Standard treatment (in 



26 



mathematical applications) employs truncating them at the optimal order growing 
with energy (or even using a Borel-type resummation if the analytic properties are 
known a priori). Instead, we adopt a cruder, but more physical definition of the 
truncated OPE series: we simply discard the OPE terms with powers suppressed 
beyond I/wIq. That is, we keep only a fixed number of powers, the first few terms in 
the OPE. Although including higher terms may further improve the approximation, 
it is superffuous for our purposes and is not adequate having in mind the reality of 
the practical applications. 

To estimate the expected amount of duality violation, we need the threshold be- 
havior of the resonances and the strength of the threshold amplitudes. In principle, 
in the heavy quark limit the least suppressed are the transitions to the P-wave states 
with amplitudes {P n \J,j\B) oc v. This yields structure functions w\, u>2 proportional 
to k 2 oc A 2 where A = m 2 — M 2 determines the energy above the threshold of the 
resonance. However, duality violation is governed by the transitions to resonances 
with Mx — Mb, for which the heavy quark limit is not applicable. We should there- 
fore assume that the threshold amplitudes do not vanish. The minimal possible 
threshold suppression comes from the phase space \k | and from the lepton tensor 
k 2 g fiu — k fl k u . Upon integrating over k 2 in Eq. (|T/|) this gives the leading threshold 
contributionQ oc A 5 . Thus, the relevant case is 7 = 5 in the notations of the Ap- 
pendix. (We also note that dimensional counting yields k = 6 for the scaling of the 
relative threshold width: ^£ ~ , A ^ 6 •) 

Using the analysis given in the Appendix, we find 

w~°' 004 ra' (40) 

where BRi ast denotes the branching ratio of the last fully open principal resonance 
in the perturbative continuüm. The latter rate is to be evaluated as yielding, upon 
averaging over the successive resonance mass gap, the perturbative decay rate: 

r M l d££-t j M 2 

V tat * S "l-^ m * , (41) 

6 T B 



AM'i 



where Ar denotes the resonance spacing. The ö(a s ) perturbative spectrum 
is explicitly known. Instead, we can use the simple dipole radiation approximation 
where the gluon radiation probability is given by |27j 



1 drpcrt „ ^M v 2 1 U2) 

rdJl4 ~ 9tt 2M x {M x -m c )' V ; 



11 As has been emphasized before, duality in the total widths depends 011 the amplitudes at k 2 = 0, 
i.e. 011 the threshold bchavior of the total (rather than differential) width. Therefore, \k\ ~ A, and 
not vmfcA. 
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with the dipole radiation coupling known to two loops and v the effective velocity 
of the charm system (recall that q 2 ~ 0): 

Ml -Ml 

Such a relation would be only approximate in the actual B X c iv decays since 
Md is not large on the scale of the radiated gluon energy. Therefore, we adopt 
an even simpler direct phenomenological bound. Since the invariant mass near the 
threshold Mx — 5GeV is significantly larger than Mp, we assume the usual linear 
in mass square spacing of the successive resonances. Next, we generally attribute 
the domain of Mx above M pert ~2.5 GeV to the perturbative excitation probability. 
The integrated rate in the numerator of Eq. ([H]) is then estimated as the total 
perturbative fraction T^ rt multiplied by the ratio of the mass square gap to the 
total length of the perturbative interval in M\: 

1 AM 2 r pcrt 

6 BRlast = Ml-MlJ^YjB) (44) 

This is clearly an overestimate, since the high-mass end of the spectrum is suppressed 
by the decrease both in the lepton phase space and in the effective QCD coupling 
at large energies. Yet this is justified since we strive to obtain an upper bound for 
the duality violation. 

For usual quark- antiquark states the gap AM 2 constitutes about lGeV 2 . As 
argued in the beginning of the section, appreciable duality violation can emerge if 
this mass gap is larger for gluonic excitations. Therefore, we take AM 2 as largeQ 
as 5 GeV 2 . 

Assembling all pieces together, we arrivé at 

-ppert 

|X| £ 0.006-^— (45) 

ppcrt 

The fraction of the perturbative semileptonic width r ^' constitutes about 10%. 
This direct experimental estimate is in agreement with the theoretical value derived 



from Eq. (j42|) . Thus, taken at face value, we arrivé at a maximal amplitude of the 
duality violation in this mechanism only at a permill level: 

Ixl £10" 3 . 

To be conservative, we can increase this estimate by a factor of three to five mak- 
ing up for the approximations we have made, most notably using the asymptotic 
expressions and neglecting the subleading in l/m& components in the amplitudes. 
We still end up with the amount safely below a half percent. Alternatively, we can 
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Clearly, such a high scale could not but affect noticeably the perturbative corrections to R T . 
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model the preasymptotic effects by assuming softer threshold suppression. The sig- 
nificant numerical suppression of \ comes from the high fifth power of the threshold 
energy gap. This power proper ly deter mines the eventual suppression of the duality 
violation; however, with high power one may have to pass more resonances before 
the asymptotic counting rule sets in to work. To account for this, we can reduce 
the effective threshold index 7 by a factor of 2. Taking 7 = 5/2 (this would, in 
a sense correspond to considering the massive lepton pair contribution which may 
dominate in the domain of intermediate energies) we would have a softer numerical 
suppression factor |% | still yielding |%| <^ 1.5 ■ 1CT 3 . 

Since the resulting numbers may seem to emerge surprisingly small, it is worth to 
address the qualitative breakdown of the suppression factor. The most obvious factor 
is the overall perturbative suppression, 0.1 to 0.15. Since the thresholds at Mx ~ 
M B ~ 5 GeV lie in the perturbative domain at least qualitatively, such a suppression 
is unavoidable. An additional factor is the fraction of this width attributable to a 
single channel. It is rather uncertain, however it emerged only 0.25, and it clearly 
cannot significantly exceed this amount. 

The most significant suppression about 1/40 comes from the nature of the duality 
oscillations, the most general fact that duality violation can only oscillate around 
the average predicted by the OPE. While the related suppression is obvious as long 
as the theoretical applicability of the OPE to the inclusive widths is not challenged, 
its magnitude may not be fully appreciated a priori. In fact it is quite significant 
and rests simply on unmodifiable mathematical constants. Let us illustrate this fact. 

Suppose we start with the most singular meaningful threshold behavior given 
by (E — E n )° 6(E — E n ). Simply subtracting the average already makes the local 
duality violation suppressed by a factor of 2 - it is described by the function (|— 2), 
< x < 1. Increasing the threshold exponent 7 by unity would amount, according to 
Eq. (^) to roughly multiplying the maximal deviation by a factor (7 + l)/27r, with 
2-7T simply being the "wave vector" of the lowest Fourier mode on the unit interval. 
Thus, already for 7 = 1 one has the numerical suppression 1/An ~ 0.1, and with 
7 = 2 it becomes l/Air 2 ~ 1/40. 

It is evident that these principal suppression factors are unavoidable as long 
as one has entered the higher resonance domain. With the numerical suppression 
for the threshold behavior ~ (E — En) 1 already amounting to more than an order 
of magnitude, it is clear that no effect exceeding a fraction of the percent can be 
identified with such a mechanism. 

As explained above, this conclusion, although not having the status of a the- 
orem, nevertheless rests on quite general arguments and relies only on very mild 
assumptions. 

A prominent effect which typically further suppresses the duality violating oscil- 
lations is the finite resonance width, which has been completely neglected. We note 
that due to suppressed threshold behavior of the semileptonic widths, the oscillations 
are already rather smooth. Therefore, the effect of the resonance mass smearing is 
far less radical than, say, in e + e~ annihilation (the case similar to 7 — > — 1 in 
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our notations) (cf. Ref. pSfl), as long as the resonances remain reasonably narrow, 
iiT n <C AM n - the damping of the oscillations is given by exp (-^jj-f-)- 

7 Comments on the literature 

Calculability of the inclusive decay widths of heavy flavor hadrons in the short- 
distance I/tuq expansion has been questioned from the time when the OPE-based 
methods were elaborated and applied to quantify the nonperturbative effects in b 
hadrons [O, [01 ITSf. The arguments in the 1990's simply contended that the OPE 



as an expansion in powers of l/mq is not applicable to the total decay widths in 
presence of confmement. However hardly any concrete evidence was offered for such 
claims nor any faults in the derivation of the OPE pointed out. 

After the validity of the OPE for decay widths had been illustrated explicitly in 
a number of nontrivial cases, the prevailing terminology has gradually changed, but 
not necessarily the substance of the criticism. It is raised by referring to "quark- 
hadron duality" as an "additional hypothesis". Yet closer scrutiny reveals that 
whenever concrete objections are formulated, the alleged effects are a violation of 
the OPE rather than of duality. The arguments thus do not differ much from 



claims in Ref. |29[ of disproving the OPE for inclusive widths based on numerical 
computations of the decay widths in the 't Hooft model. As is well known, a direct 
analytic summation of the decay widths in the 't Hooft model demonstrated full 



agreement with the calculated OPE |TÖ|, [17]]. Here we brieffy comment on more 



recent publications referred to as "supporting a possible sizeable source of errors 
related to the assumption of quark-hadron duality" |3(| and raising doubts in the 
"quark-hadron duality ansat£ J |3(J . 



7.1 "Duality violation" through 'inapplicability' of the OPE 



The first paper |JTJ considered a toy model for the semileptonic b — > c width in 
the SV limit where the energy release barely reaches the resonance domain. As 
mentioned above, the SV limit is a convenient playground for inclusive widths. 
However, to gauge the validity of duality the width obtained in the model was 
compared to the parton-level free quark width rather than to the OPE series - 
i.e., taken without any calculable nonperturbative corrections. Even the constraints 
for Mb — rrib following from the "optical" sum rule [3^] were not observed. (It is 



worth noting that in a few studied examples including the leading OPE corrections 
decreased the magnitude of the difference by more than an order.) Yet the difference 
that was found was interpreted as a violation of the OPE in the cases with 'hard- 
confined' quarks]^] It is difficult to accept such an interpretation since the OPE 
does predict such a systematic difference! 

13 One of us had discussions with the author about these matters after Ref. |?ÏJ appeared, and 
we have reasons to bclicvc that he is not fully bchind such a claim anymore. 
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Another recent paper also considered SV semileptonic decays in a toy model. 
It is more elaborate and addresses higher powers in the expansion. It claimed to 
identify a new mechanism of duality violation associated with would-be contributions 
to the inclusive width from resonances located above the kinematically allowed mass 
range in the final state. One should keep two facts in mind: Firstly, this claim 
actually refers to an explicit violation of the OPE and, as such is not relevant to 
the violation of local duality per se as it has been consistently defined. Secondly - 
contrary to how the paper is sometimes quoted - it does not claim any numerically 
significant effect. 

The basic idea behind the conclusion of Ref. is actually the same that drove 



the thrust of Ref. - that the OPE wrongfully picks up contributions associ- 
ated with transitions to states kinematically forbidden at a given b quark mass. It 
happens since the OPE allegedly is insensitive to the threshold factors 9(Mb — M n ) 
expressing, in a sense, the conservation of energy in the on-shell processes. For 
naively 6(Mb — M n ) can be dropped for any particular n at m& — > oo, or at least 
is not literally expandable in l/m& or 1/E r . In the SV approximation the partial 
width of a state P n is proportional to (Mg— M n ) 7 (generally, to an odd power 2N + 5 
of (Mb — M n ) in the model of Ref. ||33|| ). Therefore, the extra "duality-violating" 
contribution emerges always negative, and the actual width, accordingly, would have 
always exceeded the OPE expression, by an amount fading out as a power of energy 
release. 

The last simple observation shows that this would signify a violation of the OPE 
for the width rather than a manifestation of local duality violation which has to 
oscillate around zero. The total absence of such an effect is actually easily illustrated 
in a complementary way For one could apply the very same reasoning to weakly 
coupled non-confined quarks (which is how they appear in perturbation theory). The 
spectrum of "hadronic excitations" is then continuous, the 'kinematically forbidden' 
states are there at arbitrary energy release, and they even are not power-suppressed, 
only by powers of the coupling a s (for more details, see |34. 12. 27, E3|). According 



to Ref. JEJ, the OPE would completely fail - while everyone would agree it works 
perfectly in this case. Actually, here the violation of local duality would totally 
vanish since there are no thresholds and the perturbative diagrams are smooth. 



Since Ref. |33] is clearly formulated, it is not difficult to identify the subtlety 
which led to the erroneous conclusion. It lies in the fact that there was no actual 
OPE performed, but instead a simplified version was used which is only illustrative 
in elucidating physics behind the OPE results for the decay widths and its relation 
to the sum rules for the heavy quark transitions 



It can be easily shown that the 1/E r expansion adopted in Ref. |33| coincides 



term by term with the OPE if the masses of all final state excitations do not exceed 



a certain mass (Md + /^o)- However, in this case Ref. [33 would observe exact 
vanishing of duality violation at E r >fio, in full accord with our general arguments. 
We also note that the OPE in 1/E r for the transition amplitude would have had a 
finite radius of convergence given by 1//xq; at E r < /xq the OPE series simply diverges, 
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even though the formal expansion of the width itself has only a few terms. 

If the residues of the high-mass resonances do not vanish completely but, never- 
theless are exponentially suppressed, 

Pl ~ e~ cn0 , (46) 



the 'practical OPE' series still coincide with the expansion of Ref. ||33|| . However, 
the "condensates" grow exponentially, and the series for the transition amplitude is 



only asymptotic, much like in the illustrative cases considered in Refs. [y, |ÏÖ| . What 
is crucial, however, is that the possible duality violation is exponentially suppressed 
as well, like in any other Euclidean quantity. There are no specific violations of local 
duality. 

A defmite violation of the OPE was found in Ref. in the case where the 
residues p\ fade out only powerlike, 

p\ ~ (M n -M D )~ l , M n -M D ~ n a . (47) 

However, the employed expansion in this case does not coincide with the OPE start- 
ing just the power with k = l/a + 1 where the discrepancy was found. Actually, 
the expansion simply cannot be extended to this order since the corresponding mo- 
ments of the SV structure functions would diverge even in the toy model. The 
actual OPE, however still exists at this level, and yields the right result! Let us 
briefly explain this. 

In the OPE, one considers the transition amplitude A(u>) of Eq. ( |Ï2|) and expands 
it in l/(uj — E r ). For semileptonic decays one can use the representation of Eq. ( |2Ö| ) 
to get closer to T(q , q) of Ref. (f33|), and again the large-u; expansion is examined. 



This corresponds to the asymptotic expansion of T(q ,q) in l/(E r — q ). It is an 
easy exercise to check that to a first few orders in l/(E r —q ) the amplitude has the 
expansion 

T ^ = 1^— + + Tl^Sïï + T^Sl + -)• (48) 

E r -q \(E r -q ) 2 (E r -q ) (£ r -g ) 4 / 

Integrating the discontinuity of this expansion over q would exactly reproduce the 
moments through order k. (The support of ImT for the expansion above lies at 
E r — qo — 0, so one can integrate over the actual physical interval, or formally till 
infmity.) 

However, the terms starting with become different: they acquire nonana- 

lytic pieces like 

ln(g -^) 1 ln" (q -E r ) 

(q -E r y ' (q -E r y+P ' (q -E r Y+P 1 ] 

with some (generally, non- integer) powers of E r — qo or logarithms. Their disconti- 
nuity spans all the way to infmity and simply cannot be integrated over the whole 
half-axis. In the OPE we integrate it up to the proper threshold value go — "ïfc— rn c , 
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rather than to infinity as in Ref. Thus, in the actual OPE the difference 

r(B) — Tope is not given by the integral over the contour Cl along the kinemati- 
cally forbidden range of masses, as was ad hoe postulated ab initio in the paper. 
Speaking concretely, the ansatz of Eq. (|47|) 



p 2 n ~r(M n -M D y\ M n -M D ~mn p 
through dispersion relation leads to the following asymptotic term in T(qo,q 



7Ï 



y 2 (50) 

m \ f3sm(l-±)7T ( qo _ Er y-ïi+ 1 
(at physical q in the decay qo — E r < 0); at integer l — -^it becomes 

<2 r \n(q -E r ) 



v 



(3rrè (E r -q ) l ^ +1 



Moreover, the residues p n fade out powerlike only since a large momentum scal- 
ing with M n — MD can be exchanged there with the light cloud. Correspondingly, 
the related terms are not missed and do appear in the short-dist ance OPE of the 
heavy quark Green function in the form of the high-dimension operators, with the 
coefficients containing the gluon coupling. 

It is not difficult to check explicitly that the actual OPE does correctly reproduce 
the terms (and a few more) in the exact width T(B). Instead of presenting the 
straight forward derivation, we give here a simple heuristic argument which is more 
than convincing. The point is that the amplitude given by the leading term(s) in 
expansion (|49|) is very similar to the effect from usual gluon perturbative corrections 
with nonconfined quarks and gluons, however with the 'coupling' suppressed by a 
corresponding power of the energy. As explained above, nobody expects the OPE 
to fail and local duality to be violated in such a "non-resonant" situation. This 
simply means that going through the formal OPE machinery, viz. integrating the 
imaginary part of the asymptotics of the transition amplitude is bound to yield 
exactly the averaged contribution of the high hadronic thresholds. The average 
difference between Tope and T(B) would vanish. 

To summarize : the observed violation of duality of Ref. has its root simply 



in the fact that the actual OPE was replaced by a procedure differing just in terms 
found as an alleged violation of duality. 

We note that the most recent paper by the same authors |2(J appeared while 
this Vademecum was in writing, dedicated to the question of duality. It does not 
anymore claim establishing duality violation at the considered level. We view this 
as an indication that the authors do not stand behind the mechanism of duality 
violation of Ref. p3| |. Moreover, their analysis clearly supersedes the discussions of 



Isgur jïï] , and that chapter can be considered closed. Yet we note a question raised 
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in the footnote about the validity of some conclusions of Ref. \TÜ\. According to the 
above illustrations, such qualifications should be regarded as groundless. 



A general note can be made here regarding the attempts to challenge the ap- 
plicability of the OPE. There are two basic ingredients involved in the 'practical 
OPE'. One is constructing the large-a> expansion of the transition amplitudes like 
A(uj) in Eq. (|T|) in terms of the local heavy quark operators. The second step is 
relating this asymptotic expansion to the actual width. It is understandable that 
the challenge is motivated by the (sometimes indeed not very intuitive) fact that the 
final expression for the total width is given by the expectation values of the local 
heavy quark operators. However, it turns out that all the concrete objections in the 
literature or those hanging around as folklore attempted de facto to undermine the 
second step, which is a purely mathematical procedure and is indisputable. 

It is the first step which is far less trivial, both physically and technically. Namely, 
here we relate the high-energy asymptotics of the transition amplitude in QCD to 
the short-distance expansion of the heavy quark Green function in the external 
field averaged over the nonperturbative QCD ensemble of quarks and gluons (recall 
Fig. 1). And, curiously enough, just this least trivial part has not been challenged 
in the publications attempting to disprove the theory of the 1 / mg expansion. 

7.2 Local duality proper 

We do not have to say much regarding the literature in this respect, for - to 
our knowledge - motivated estimates of the significant effects in total semileptonic 
widths have never been presented.Q The realistic estimates we mentioned all yield 
tiny duality violation. The possible rationale has been illustrated in Sect. 6.2. 

8 Quark- hadron duality and extracting \V c b\ 

One of the central phenomenological applications of the heavy quark expansion in 
beauty decays is extracting the value of the CKM parameter V c b. Two quantities 
are believed to be best suited for this purpose yielding the highest precision with 

1 The situation at times contains an ironie twist when habitual critics of local duality in the 
intcgrated semileptonic widths present quite precise relations of thcir own on partial spectra in 
semileptonic or radiative B decays - with such concerns becoming muted. For instance, it has been 
put forward (see, e.g. |35|) that a measurement of part of the photon spectrum in B — *■ X s + 7 
enables us to predict the B — > X u iv rate restricted to a small and most vulnerable slice near the 
end point with a few percent accuracy. One has to keep in mind here that such a rclation ]36| 
rests solely on the OPE and actually holds only to the leading order in l/mg. It is not otherwise 
supportcd by heavy quark symmetry or any other independent argument. 

Similarly, reservations regarding the accuracy of duality for total semileptonic widths can be 
found in Ref. |57j, while estimating in Ref. |38| the uncertainty in the b —> u rate over the very 
limitcd (and marginally hard) domain of maximal q 2 to be only a few percent, including the effects 
of local duality. 
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limited model dependence. One is the total semileptonic width r s i(6 — > c); another 
is the exclusive rate for T(B — > D*£u) extrapolated to zero recoil, where the B — > D* 
formfactor F(0) is fixed if both and m c are sufficiently heavy.0 

The experimental advantages are well known: the total semileptonic b — > c rate 
is by now one of the best measured quantities in B physics. Measurement of the 
B — > D*£u rate is far more involved, but due to the pioneering work by ARGUS 
and CLEO, and also with LEP and recent BELLE data has been determined with 
decent accuracy. Yet the necessary extrapolation to the point of zero recoil still 
introducés sizeable uncertainties.0 The recent significant jump in the CLEO value 
of the zero-recoil decay amplitude is an illustration of these problems. 

Here we address the purely theoretical aspects of the comparison of the two 
methods to sort out facts from fiction concerning these approaches (a brief review of 
the historical misconceptions can be found in Ref. P6|). In particular, the question of 
'quark-hadron duality' has been advertised as an additional assumption inherent 
and specific to T S \(B). 

Once again, we have to stress that the OPE expresses the width as a power series 
in the inverse energy release 

_M = O0 + a2 _ + a3 _ + .... (51) 

All such series are asymptotic and do not define the sum with unlimited accuracy 
at finite masses. As pointed out by Shifman 0, it is this feature of the expansions 
that underlies violations of local duality. Clearly the same is true for the B — > D* 
formfactor at zero recoil: 

F(0) = b + b 2 \ + b z \ + ... , (52) 

m ó c 

where we have dropped the numerically less significant 1 / mj corrections. This series 



is likewise asymptotic - the expansion Eq. fl52 ) does not define completely -F(O) as a 
function of quark masses. The difference in respect to general quark-hadron duality 
between the two cases, therefore lies in the details of the expansions potentially 
affecting the quantitative behavior. 

The two series of Eqs. fl5*T|) and (|5*2"| ) superficially look quite similar; most notably, 
neither contains a correction linear in I/tuq. Yet there are substantial differences. 
The expansion for T S \(B), Eq. (|5"ID , is given by B meson expectation values of local 
heavy quark operators. From a theoretical viewpoint, they are of a rather universal 
nature, and their size can be related to other short- dist ance observables. In the 
series for -F(O), Eq.(^2[), on the other hand, none of the power corrections is related 



to local operators. Instead they are expressed by typical long-distance correlators 



15 Both methods were discussed and even applied already in the 80's. Their actual value has 
been revealed with developing the theory quantifying the preasymptotic corrections [f26| . 

16 Thc claims of a rcliable theoretical prediction of the shape of the B — > D* formfactors quoted, 
c.g. in Ref. (ÏÖ), cssentially overstate their alleged model-independence. 
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shaped by dynamics of momentum scales ~ ö(Aqcd)- Even the leading power terms 
are basically unknown. It should be noted that the emergence of any such term in 
the series for T S \(B) would constitute a breakdown of the OPE and invalidate the 
whole theory of the inclusive widths. 

As stated above, duality violation is related to the effects of sufficiently high 
orders which for asymptotic expansions may not fade out quickly enough; including 
higher orders will not improve the accuracy of the result. While this represents a 
legitimate question for Eq.(^I|), the question itself is somewhat ambiguous for -F(O) 
in Eq.([52|), since there already the first terms are unknown in practice. 

Nevertheless let us ignore this practical problem and just assume that all terms 
in Eq.(|52|) can be determined with unlimited accuracy. Then we can address the 
issue of duality for -F(O) and compare it with the situation in T S \(B). 

According to conventional wisdom the series ([52]) approximates the exact ampli- 
tude with exponential accuracy for nearly Euclidean quantities like F(0): 

5 D F(0) = F(0) - F(0) 1/mcexpansion ex e- mc/M " , (53) 

with Mh representing the characteristic hadronic scale (for clarity we keep only the 
leading corrections related to m c ; i.e., we put — > oo while keeping m c finite). 

Such an exponential suppression of the deviation is not necessarily the case 
for a Minkowskian OPE observable like T S \(B). Even though the appropriately 
averaged widths coincide with their expansion exponentially in E r , the point-to- 
point oscillating difference can in principle be only power suppressed: 

£12-1= <**)>=„. (54, 

1 OPE &r 

Thus, the assumption of local quark-hadron duality is the approximation that one 
can discard, without averaging but at a given large value of the b quark mass and 
E r the power-suppressed r.h.s. of Eq. (Q) like the "usual" exponential terms in the 
r.h.s. of Eq. (0). 

Because of the different functional behavior we naturally expect the OPE to ap- 
ply better numerically in the Euclidean than in the Minkowskian domain: there may 
be the oscillating component in inclusive probabilities. However, this general obser- 
vation obviously requires a crucial qualification: the observables and, in particular, 
their "hardness" in energy scale must be commensurate. 

In practice, however, the situation is rather different in -F(O) vs. T S \(B). The 
latter is much better for the OPE than a generic inclusive probability, say than R T 
(not to mention R e + e - or T n \(B)). Even more importantly, the large mass expansion 
parameter is far better in r s i(-B) (EV ~3.5 GeV) than in -F(O) (m c ~ 1.25 GeV). This 
underlies the actual hierarchy. While at large enough m c the exclusive F(0) would 
evidently enjoy a smaller duality uncertainty, in practice T s i(B) is much more 
robust in this respect. 

The exponential suppression ([x|) in the Euclidean case will actually set in only 
for masses that sufficiently exceed the scale of the strong interactions. Duality 
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violation in the semileptonic width is numerically too strongly suppressed in this 
case to be relevant in practice regardless of being exponential or only powerlike. We 
illustrate this point in numbers. 

Let us recall that the leading long-distance term in -F(O), Eq. fl5"2"|) constitutes 
an about —7% effect [p4fl ; i.e., the correction in the decay rate is —15% (cf. —5% 



in T s i(B)). What can be the magnitude of the exponential terms Eq. (^) which 
are not captured in the straightforward l/mg expansion? With m c = 1.25 GeV 
the size of ó~dF(0) would constitute dozens percent even for the moderate value 
M = 600 to 800 MeV. However, we can argue that the actual high momentum scale 
set up by quark masses is here 2m c rather than m c itself.Q Then one has 

f 1.5% /ihadr = 0.6 GeV 
5 D F(0) oc e~ 2mc/ ^ ~ l 4% /i ha dr = 0.8 GeV (55) 

{ 8% /i hadr = 1 GeV 

While this simple estimate is only tentative, it clearly shows that the l/m c expansion 
has intrinsic uncertainties here at least at a percent level even in the most optimistic 
scenario. This is in accord with the realistic estimates of the magnitude of the second 
and third order effects, ö 1/m 2F(0) ~ -7% U, \S 1/m3 F(0)\ > 2% 0, |2|.0 

The estimates of SdF(0) can be confronted with the local duality violation in 
the semileptonic b — > c decay width. At the total energy amounting to the charm 
mass 1.25 GeV it may naturally be large. However, the key fact - the significant 
energy release E r ~ 3.5 GeV ^> m c makes things completely different. All estimates 
given in the preceding sections as well as in other dedicated analyses, yielded the 
magnitude much below a percent level. 

It should be noted that moderate uncertainties in F(0) are obtained only in 
optimistic scenarios where /ihadr does not exceed 0.7 GeV. If that indeed were to 
reflect the true scale of nonperturbative dynamics, then violations of local duality 
would simply be invisible in T S \(B). 

As first explained in Ref. |26| , the scenario with /ihadr — 1 GeV or even somewhat 



larger is not ruled out so far phenomenologically, even though there is no compelling 
direct evidence in its favor. To find possible significant duality violation in r s i, we 
allow the relevant /ihadr to be as large as 2 to 2.5 GeV. The quark-hadron duality 
'violation' SdF(0) does not need to blow up literally in this case, of course: even in 
the absence of a real exponential (in m c ) factor it may be suppressed by numerical 
factors of order 1. For instance, at m c <^ Aqcd the corresponding zero recoil form- 
factor for B — > K* is expected to be 0.5 to 0.6. Even in the worst scenario one still 
might have F(0), say around 0.8. The point is that in such a case the estimate of 

17 This is related to the position of the other singularities in the transition amplitude. This 
feature is explicit in the concrete lowest order calculations. 

18 The idca of insignificant higher order contributions to F(Q) has its roots in the original es- 
timates of this formfactor by Neubcrt being very close to unity: Si/ m 2F(0) = (—2 ± 1)% jtÖ|. 
Such scale of the leading nonperturbative effects would indeed suggest negligible higher-ordcr and 
exponential terms. 
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F(0) in the I/mg expansion would be not much more sensible than trying to obtain 
Fb~>k* i n a l/m s expansion. 

This discussion shows that contrary to popular lore total semileptonic widths 
suffer considerably less from duality violations than the zero-recoil B — > D* rate 
for the actual values of m& and m c . The violation of local duality inherent to decay 
widths, is strongly suppressed by the peculiarities of the total semileptonic width, 
and the numerical hierarchy is completely reversed in favor of Y S \{B') by the signif- 
icant energy release in b decays, which is much higher than the charm quark mass. 
Some damping exponents are still larger that other "slow powers" . 

9 Conclusion and Outlook 

In this paper we have critically reviewed the salient aspects of treating nonper- 
turbative QCD relevant to quark-hadron duality, and in particular to violation of 
local duality as it applies to the total semileptonic widths of heavy navors. A new 
underlying mechanism was addressed which has been suggested previously to be 
potentially significant, but not analyzed before. Still we have found in Sect. 6 that 
conservatively it cannot exceed a fraction of percent level in T S \(B) with realistic 
estimates yielding even much smaller numbers. Thus, the corresponding uncertainty 
is totally negligible in practical extraction of Va- 

It should be clearly realized that this conclusion applies only to the totally in- 
tegrated semileptonic widths. Most of the suppression factors we discussed do not 
apply even qualitatively to partially integrated widths, in particular restricted to 
the limited regions of the lepton phase space. This reservation is supported by the 
explicit model estimates. For instance, in the 't Hooft model the local duality vio- 
lation in the differential distributions, in particular, in b — > s + 7 or b — > u Iv was 
found to be very significant numerically [E2^] for practically relevant kinematics. 

Typically, we expect duality violation to affect T S \(B) only at a permill level. 
However, this does not mean we can compute in practice T S \(B — > X c ) today with 
such an accuracy. For we specifically addressed the violation of local duality proper 
here, and considered only the related uncertainties. The theoretical progress in the 
recent years allowed us to formulate the question consistently; this is far from trivial 
and has often been missed in the literature. In simple terms, this uncertainty tells 
us how accurately we can determine the width if as many as necessary terms in the 
(practical) OPE are known with sufficiënt precision. 

From the present and preceding theoretical studies of quark-hadron duality and 
its violation we can draw a final practical conclusion. The accuracy with which 
\V c b\ can be extracted from r s \(B) (and \V u b\ from T S \(B —>■ X u )) is actually limited 
only by the precision with which we know the first few terms in the OPE. Duality 
violations are negligible in practice. Even approximate knowledge of the magnitude 
of higher order nonperturbative contributions could be of help in a complementary 
respect: to determine the mass scale /ihadr which provides the effective yardstick for 
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nonperturbative dynamics in heavy quarks, /ihadr < 1 GeV would imply a "favorable 
scenario" where one can hope to apply heavy quark expansion to quantify deviations 
from the heavy quark symmetry in charm. If it were to exceed 1 GeV, one would 
probably have to abandon such methods for model-independent determination of 
underlying weak decay parameters in general, and rely only on the more robust 
cases of the inclusive observables in beauty decays. Even such a scenario would not 
elevate local duality violation in T S \(B) to practically significant levels. Yet it would 
affect extracting the fundamental input parameters like quark masses and the first 
nonperturbative parameters fj%, Hc, Pr, from the data. As we point out below, there 
is a direct experimental way to infer the scale of /ihadr i n B decays. 

All considered mechanisms yield quite small violations of local duality in T s i(B — > 
X c ) when compared to the impact on R T from r decays, which has become the 
canonical yardstick for judging the OPE. This had actually been anticipated. The 
qualitative arguments can be found in Ref. [ Ï2| ; now we can make the comparison 



in a more elaborated manner illustrating this in quite general terms. 

• The one argument given in the literature to substantiate smallness of local 
duality violation in R T invokes a strong threshold suppression, oc A 2 : the hadronic 

states X with Mx = M T — A contribute to R T with the weight (]jf-) • Yet the 
threshold suppression is much strong er in T s \ where it is the fifth power, A 5 . 

• When decreasing the energy scale R T blows up and duality is 100% violated. 
T S \(B— >X C ), on the contrary, exhibits very accurate local duality even with shrinking 
energy release; this is ensured by the SV limit. 

• Local duality violation in T S \(B— >X C ) comes from highly excited charm states 
with Mx c ~ Mb\ they are thus driven by perturbative gluon exchanges. In their 
absence the yield in the relevant domain - denoted by P in Fig. 3 - would practically 
vanish. The relative scale of such an effect starts at the a s /n <^ 15% level, even 
before any other suppression factors are considered. In contrast, duality in R T or 
R(s) affects already the "valence" quark contribution which is roughly a constant 
N c at arbitrary energy. 

• The energy scale itself is obviously much higher in T S \(B — > X c ) than in r 
decays. 

With all these effects acting in the same direction, one ends up with a very small 
local duality violation in r s i(B— >X C ) even allowing for it to be as large as about 5% 
in R T . In view of this comparison, discussing possible significant duality violations in 
the semileptonic beauty widths and not allowing for appreciable effects in r decays 
does not seem to constitute a consistent application of QCD. 

The suspicions towards large effects of local duality violations in the decay widths 



19 We remark here that - contrary to common lore - increasing this power in rcality deterioratcs 
duality: for a large power n R T is saturated at s <J M^/n, and eventually it simply ceases to be a 
short-distance quantity. Suppression of local duality violation with increase of n would only apply 
in the academie case where scales with n so that M^/n remains much larger than ^^adr- ^ n 
practice half of R T comes from s in Eq. (fj]) below 1 GeV. 
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were traditionally fed by reported problems in nonleptonic widths, possibly showing 
up in the size of BR s i(5) and of the t^/tbo ratio [p^j . However, the maximal mag- 
nitude of local duality violations in nonleptonic decays is not limited, in principle, 
by any of the features peculiar to semileptonic widths as discussed in Sects. 4 and 
6, which tend to reduce duality violation there by orders of magnitude. 

Another potential problem discussed is the absolute semileptonic decay rate of 
D mesons: about a third of the actual width may seem to be unaccounted for in the 
l/m c expansion; this deficit could be interpreted as due to duality violation. How- 



ever, as pointed out in Ref. [f|l|, the excess of the observed decay rate can naturally 
be explained by usually discarded nonvalence four-quark expectation values; they 
would still yield a marginally noticeable correction in total B decay widths. This 



conjecture was reiterated in Ref. |22| . In any case, it is natural to expect sizeable, 
or even large duality violation in charm while having it very small in beauty. 

A rather direct experimental check of local quark-hadron duality and its onset 
in a broad enough range of energies would be a measurement of the decay rate 
distribution over the invariant final state hadron mass Mx (or, equivalently, of a 
combination of the structure functions w\ and W2) in the B — > X c tv decays in the 
mass region above 3 GeV.0 This would provide unique and detailed information 
on QCD in the transition from the nonperturbative to the perturbative regime, at 
the level of radiative effects directly in Minkowski space. The primary goal here 
would be to analyze the onset of the perturbative regime, domain P in Fig. 3, by 
comparing the actual hadronic yield with the perturbative result. (The structure 
function W3 becomes accessible if the electron energy is additionally measured.) 

The total semileptonic width corresponds to integrating these structure functions 
over the whole perturbative domain and also dilutes them by adding the dominant 
contribution from the resonance region R. This makes duality violation invisible in 
T s i(i?). However, local duality violations are expected to be evident in the differential 
distribution, and it should be checked that the structure functions oscillate around 
their perturbative expressions, in particular in the lower part of the interval in Mx 
- but coincide with the latter upon averaging. As a byproduct of such an analysis 
one can obtain an independent direct measurement of the effective QCD coupling a s 
at low energies, and evaluate the magnitude of the higher-order local heavy quark 
operators entering heavy quark expansion. 

A similar - although, probably, more remote experimentally - possibility is to 
study in detail the similar Mx distribution in b — * u transitions - much like it was 



proposed for model-insensitive extraction of \V u t\ [0. Here even a more detailed 
information on the duality onset is accessible in principle via genuine doublé or 
triple distributions. In practice, though this may be obscured by the necessity 
to restrict kinematics to exclude B — > X c iv decays and by intervention of the 
Standard nonperturbative OPE corrections (in particular, due to Fermi motion) up 

20 Thc resonance region below it essentially determines the OPE parameters A(/x), /x 2 , /Uq, p'jj 
etc. 
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to M Xu = Ito 1.5 GeV. 

In general, we think that all practical methods of extracting V uh suggested so 
far would bear a more or less significant element of model dependence until local 
duality is carefully studied in B — > X c lv along the lines suggested above. 

While the central role of the OPE is becoming widely accepted now, certain 
subtle, yet significant features of its implementation are often not fully appreciated 
despite extensive reviews and successful applications in the literature. It is often 
overlooked that only a careful treatment of the OPE yields accurate and defendable 
results: 

X One should use properly defined cutoff-dependent strong interaction param- 
eters: quark masses m& )C (/i), kinetic fJ^(fj) and chromomagnetic /^(/-O expectation 
values, the Darwin term p 3 D (fi) etc. Ill-defined and indefinite pole masses, Ai of 
HQET etc. should be avoided. 

X In this approach there are no significant uncertainties in the radiative cor- 
rections or those associated with the value of the b quark mass. This has been 
substantiated by explicit computations of the full ö(a 2 s ) corrections [jï3 |. 

X There is limited direct impact from the uncertainty in /i 2 and /Xq, as well as 
from the higher-order terms. 

The analysis shows that the current limiting factor in determining V c b is the 
precise value of m c . At present the most accurate way is relating m c to mj 



using the expansion of hadron masses [15]: 




M B + 3M B * M D + 3M D * f4 f 1 1 

mh—m r = + 

4 4 2 \m c mj 4 

(56) 

This is the only place in evaluating the semileptonic width where we rely on the 
expansion in l/m c and thus are sensitive to the values of the higher condensates. 
Theoretically, the nonrelativistic expansion of mass is known to behave the best 
compared to, say, the expansion of wavefunctions and their overlaps determining, 
for instance the B — > D* formfactor. There exist a number of independent de- 
terminations of the short-distance charm mass m c , see, e.g. Refs. [[PJ for m c (m c ). 
Translated to the running masses used in the heavy quark expansion, they yield very 
similar numbers using just the central values of m&, fj% and the expected moderate 
values of the higher-order condensates, without even invoking the error bars. 

Recently a direct determination of the heavy quark mass and the kinetic expecta- 
tion value has been reported by CLEO |45| . Unfortunately, the preliminary version 



gives the fitted parameters A corresponding to the pole mass, and — Ai which makes 
the results generally unusable. However, using the quoted values of the hadronic 
moments themselves and the estimate of the Darwin term p z D = 0.12 GeV 3 |Ï6 one 
obtains central values for mb(fi) and l^(/J>) in a good agreement with the theoretical 
expectations [||: A(l GeV) ~ 670 MeV, /4(lGeV) ~ 0.42 GeV 2 . We note paren- 
thetically that such an analysis, as a matter of fact, is most sensitive to the value 
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of m c (more precisely, to m^ — m c ). In this respect, it is largely suited to determine 
the possible size of the higher-order terms in the mass expansion Eq. (|56|). 

To conclude: future improvements in the theoretical accuracy in V c b are expected 
to occur through better control over the higher-order terms in the meson mass 
expansion Eq. (|56|) and/or an independent precise determination of the low-scale 
running charm mass. 



10 Epilogue 

As a final note let us express how we view as quite paradoxical the way in which 
the discussion of quark-hadron duality and its limitations in B decays has - and has 
not - taken place. 

On one hand there is the exclusive transition B — > D* Iv: its formfactors Fjj* 
cannot be expressed through an expansion of local operators even at zero recoil and 
the leading power corrections ~ 0(1/ m 2 c ) are not fully known; they are estimated 
relying on natural assumptions f24] , ^6) which, however become vulnerable when 
descending to the level of a few percent error bars. Nevertheless statements of an 
unrealistically small theoretical uncertainty are typically accepted as gospel without 
much reflection; it is forgotten how much the central value stated for -F(O) has 
changed over the years. 

The psychological background behind the apparent tolerance is quite under- 
standable: Already the leading power corrections oc l/m? c in F(0) are unknown and 
the estimates of the small-uncertainty variety come from ad hoe model assumptions. 
There is no much to criticize here beyond that. 

Yet for the fully inclusive semileptonic width T S \(B), for which the consistent 
OPE has been given, suggestions of uncontrollable theoretical errors are readily 
picked up - despite the fact that dedicated theoretical analyses have given no re- 
producible sign of such effects and that the theoretical predictions have not changed 
in any significant way over the years. Part of the reason might well be that even 
the central result - the absence of ö(l/mg) corrections in the width |T3|, [TJ| - is 



highly nontrivial and becomes intuitive only within the proper approach. The OPE 
itself, while well developed, remains a conceptually nontrivial theoretical technology 
employing a number of basic tools. This provokes, in our opinion, much critical 
attention to the total widths and attempts to challenge the heavy quark expansion 
for T s i(B— >X C ) at a more pedestrian level. 
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11 Appendix 

The analysis briefly described in this Appendix is quite general in nature. There- 
fore, we do not specify explicitly the type of observable we discuss. Likewise, we 
generically denote the large energy scale parameter as E (it is mg or E r in the heavy 
quark decays; likewise, it can actually be ttiq or s = E 2 ); however, to keep in the 
context with our discussion we will refer to the short-distance observable as Tb- 

Under our convention explained in Sect. 6.2, we discard the higher-order terms 
in the OPE, and 

fr oc -^f{E-nA) (57) 

and Tope/Tc) = b + b 1 /E + ... + b k /E k (the coefficients b k can logarithmically depend 
on E). Likewise we can discard all effects in Tb which fade out faster than 1/E k , 
since they are subdominant in x(-E). The requirements stated in Sect. 6.2 are then 
sufficiënt to determine the exact asymptotic form of the local duality violation, i.e. 
the function f(E), including its overall normalization. 

Indeed, let us consider (Tb — ropE)/To on the interval between the opening of 
the two successive thresholds. The contribution of the last open channel in Tb/^o 
is given by 

° : (E-nAye(E-nA) (l + O (^] 



E kK ' v '\ \EJ. 

We do not need to know explicitly either contributions of the other, lower mass 
channels (they are expandable in 1/E), or Tope - it suffices to know that their 
contribution to Pope/Po is given by a few powers of 1/E. Expanding the difference 
of these terms around E n = nA we get a general polynomial 

4ü (ao + a^E-En) + a 2 (E-E n ) 2 + ... + a k {E-E n ) k ) . (58) 



E k 

The duality violation is then given, up to a factor c/E k , by the difference 

9l (^) = (Ai?) 7 - oo - a,AE - a 2 (AE) 2 - ... - a k (AE) k (59) 



A£ = AFrac(£) , X {E) = c^i 

(Frac and Int denote the fractional and integer parts, respectively). 
The function g 7 is constrained by the requirements that 

• its average vanishes 

• it is continuous, as well as its derivatives up to the order l = Int (7). Higher 
derivatives have discontinuity at the threshold values E = nA, determined by the 
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threshold behavior (Ai?) 7 . 

These impose 1+2 linear constraints on the coefficients a^. Therefore, the consistency 
demands k>l + l. In the cases we consider one actually has k = l + 1. This is not 
always the case; the thresholds can appear at a level suppressed by a higher power k 
than the threshold exponent 7. However, it can be argued that in this case the last 
few (viz., Int(/c — 7) — 1) terms in the polynomial are absent, since it emerges from 
expanding the terms power suppressed at least as ^- Int ( fc -7-i\ This observation 
is quite general and is related to a certain property of the OPE which is similar 
to global duality |Ï5 |. We do not pursue this aspect here, and simply make use of 
the fact that k = l + 1 in our applications. Then the whole function describing the 
asymptotic violation of local duality to the leading order in 1/E is fixed in terms of 
the large-i? threshold behavior. 

A few relevant examples are given below. First, however, we note that the 
branching fraction of the last open channel is c ^^? 7 , and it is maximal just at the 
threshold for the next resonance. Therefore, it is advantageous to normalize the 
amount of duality violation to this maximal last resonance fraction BRi ast : 

X (E) ( AE S 



X(E) 



BR 



last 



#7 



V A 



(60) 



Likewise, the maximal magnitude of such a normalized amount of duality violation 
|Xo| is simply a universal number depending only on the power of the threshold 
suppression. 

Now, for the most singular threshold behavior possible in 3 + 1 dimension A2 
we have 



f \ 1 1 
9{x) = £ 2 - - - x 

o 



< x < 1, 



|Xo| 



1 

6 



(61) 



For higher power s we get 



9{x) 
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g(x) 
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IXol 




0.0282 


at 


7 = 


3 
2 
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IXol 
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at 


7 = 
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|Xo| 
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(62) 



The case 7 = 2 is relevant for the r width. It was obtained in Ref. [ÏD|| by explic- 
itly constructing the OPE. In lower space dimensions a more singular behavior is 
possible. The case 7 = was considered in Refs. DJ for the 't Hooft model. 
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Figure 4: Functions g 7 (x) for integer 1 to 5 (a) and half-integer \ to | (b) threshold index 
7. Curves with smaller magnitude corrcspond to larger 7. Lower plots show magnificd g 1 {x) for 
7 = 3, 4, 5 and for 7 = |, |, |, respectively. 



The corresponding functions are plotted in Fig. 4. The maximal magnitude of 
the duality violation decreases fast with smoothening the threshold behavior, which 
is transparent. Increasing 7 by one requires vanishing of the average of one more 
derivative of the function. This amounts to integrating g(x) and subtracting its 
average: 

r cv. ^ 1 

(63) 



g 1+1 {x) = (7 + 1) £g y (t)dt + J tg y (t)dt 



its maximal deviation from zero drops fast. It is remarkable that the g 7 (x) for 7 > 1 
is very close (up to a factor) to cos (2ttx + 7 ); the origin is clear from the above 
recurrent relation. The latter simply states that the n-th Fourier mode of g 1 {x) gets 
multiplied by 2(7 + l)/(27m); the higher modes then die out quickly and only the 
leading one with n — ±1 survives. It is worth noting that the normalized deviations 
|Xo| do not depend even on the particular mass pattern - the resonances can be 
equally spaced in mass or mass squared, which is an added convenience. 

It should be clear that the above evaluation of the duality violation is asymptotic. 
At intermediate energies subleading in 1/E effects can be noticeable; likewise the 
asymptotic behavior of the resonance masses is distorted, which modifies the exact 
numbers. 
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